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ABSTRACT 


Water in the soil is one of the key environmental parameters. With increasing 
attempts to understand and to manage the environment, the demand for on accurate, 
timely, and locaHon*^pecific soil water content information can be expected to 
grow. Two types of techniques appear potentially capable of providing such infoP" 
motion, namely water balance and remote sensing methods. In this study, the 
attention was focused on the microwave remote sensing method, and in particular on 
bare soil temperature and moishire regimes and their impact on the microwave technique. 
A layered water balance model was developed for determining soil water contents and 
their changes in the upper zone (top 30 cm of soil), while soil moisture regime charac” 
teristics at greater depths and those near the surface during the diurnal cycle were 
studied using experimental measurements. Soil temperature and its variations due to 
several parometers were Investigated by means of a simulation model. Using the two 
models, moisture ond temperature profiles of a hypothetical soil located at mid-latitudes 
were generoted, analyzed, and subsequently used for computing microwove soil parameters 
at three frequencies (1,4 GHz, 4,0 GHz, 10,0 GHz) for a dear-sky summer day. 

The results suggest that (i) soil moisture in the top 30 cm can be predicted on a daily 
basis for 1-cm depth increments within approximately 0,02 cm /cm”^, (ii) soil temper- 
ature presents no problem in microwave remote sensing of soil moisture provided that 
surface soil temperature can be measured with infrared radiometers from the same plat- 
form, and (iii) surface layer domination exists, i.e., remotely measured microwave 
response of a bare soil is determined primarily by the moisture at and near the sur- 
face unless appreciable air-soil discontinuities (in terms of wavelength) exist within 
the soil , In conclusion, an algorithm is proposed which combines the water balance 
and microwave methods to monitor profile soil water content changes over large areas. 
Finally, a summary of passive and active microwave measurements of soil moisture is 
presented . 
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CHAPTER 1 . 


INTRODUCTION AND PROBLEM STATEMENT 
1.1 INTRODUCTION 


Soil water, due to its presence near the earth's surface, plays a fundamental 
role in mass and energy exchange processes at the lithosphe e/atmosphere interface 
and In the ecosystem's functioning. Its importance has long been recognized and is 
presently topical ' in such diverse scientific disciplines as hydrology, ecology, plant 
physiology, meteorology, climatology, geomorphology, hydrogeology, pedology, 
etc. Hydrology and ogricuiture are concerned primarily with those features of soil 
water which affect man. The physical geographer's interest in soil water stems from 
the important function surface waters have in differentially shaping the earth's 
topography, the latter occurring mainly due to spatial variations of water availability. 

In a recently proposed unifying concept of physical geography (Carter et a| ., 1972), 
soil water was considered a basic component of the interface environment. Because 
of its indispensability for survival and a prosperous livelihood, soil water also affects 
various phenomena in the domain of cultural geography. 

A necessary requirement in dealing with many aspects of water in the soil is 
some means of measuring the amounts present over large areas. The numerous methods 
developed for estimating soil water content may be divided into two groups based on how 
the moisture value is obtained. 

1.1.1 Direct Methods 

In the direct approach, water content of a soil segment at a given time is 
determined by directly examining the segment at that time. Thus discrete moisture 
values are provided, although in some cases continuous recording is possible. The follow" 
ing direct methods are used most often. 

Gravimetric Method 

The segment is physically removed, and its water content is determined by 
separating the sample into dry soil material and water components, respectively. The 
method is destructive, t.e., the soil segment is destroyed during this procedure and thus 
repeated determinations are not possible. 



Res15^ance Mehhod 

The resisl'ance to electric current flow decreases with increasing moisture content 
along the path of the current. This decrease can be measured and through calibration 
related to water contents. Blocks made of porous materials such as gypsum or nylon 
(Baver, 1956) are inserted into the soil segment to be measured. The method is non~ 
destructive . 

Tensiometric Method 

Due to the forces of cohesion and adhesion existing between water molecules 
or between water molecules and soil particles, a force inversely related to the amount 
of water present must be exerted to remove the water out of the soil segment. The 
magnitude of this force (tension) can be determined by measuring the force applied by 
the soil on a pool of free water across a permeable membrane. The measured force can 
then be related to soil water contents through calibration curves for a given soil. This 
approach yields repeated measurements of a soil segment once the tensiometer is installed; 
it is not opplicoble, however, for tensions higher than about 0.85 atmosphere {Taylor 
and Ashcroft, 1972). 

Psychrometric Method 

Relative vapor pressure of soil water depends on the difference in water potential 
that results from changes in matric, osmotic, and pressure potentials. Although the 
vapor pressure range is relatively small for water contents in the plant growth range 
(Taylor and Ashcroft, 1972), accurate measurements can be made provided that the 
ombient temperature effect on the vapor pressure is taken into account. For example, 
Rawlins and Dalton (1967) developed a field psychrometer that can measure soil water 
potentials with a precision of +50 joules/l^ in soils where ambient temperature changes 
less than 4®C , 

Scattering Methods 

When high energy neutrons collide with atomic particles of approximately the 
some size, they are rapidly slowed down by successive collisions until they attain the same 
velocity as particles with which they collide. Since the hydrogen atom is most effective in 
slowing down fast neutrons, the number of slow neutrons will be proportional to the amount of 
water present. The slow neutron count may be detained by meons of a detector inserted 
into the soil segment, and water .content within the segment con be determined from a 
calibrotion curve. This approach is called the neutron method (Bell, 1973). 
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Another way of estimating woter content nondestructive I y is by measuring 
attenuation of gamma rays as they pass through the soil segment. This method is most 
useful for laboratory work since it requires that the soiree and detector be placed on 
opposite sides of the segment (Gardner and Calissendorff/ 1967). However, the technique 
has also been applied in the field (Ryhiner and Pankow, 1969; Reginato and Jackson, 
1971). 

The direct methods supply relatively accurate data about the soil in close 
proximity to the sensor, or in the ease of the gravimetric method, about the sample 
taken; however, their applications are limited by the cost of instrument Installation and 
operation . 

1*1.2 Indirect Methods 

Since they do not require physical contact with the soil, indirect methods 
represent the only realistic means of determining soil water content over large areas. 

There are two basically different indirect methods, namely water balance and remote 
sensing. The estimates by remote sensing are discrete in time while the water balance 
approach can in principle be used on o continuous basis* 

Water Balance Method 

Rother than by direct measurement, water content is determined os a residual 
after all inputs (precipitation, irrigation) and outputs (runoff, drainage below the 
terminal depth, evaporation, transpiration) have been accounted for. While most 
techniques in this category yield water content for the entire soil column, some models 
(e.g.. Baler and Robertson, 1966) subdivide the column into several zones. Water 
balance models have proven useful and reasonably accurote (Baier, 1967). More 
successful applications are hindered by the necessity for areal extrapolation of point 
measurements which is intrinsic to the method. 

Remote Sensing Methods 

Remote sensing of soil moisture is based on the change in electromagnetic 
properties of the soil resulting from an addition of water. The changes and their 
potential usefulness for soil moisture estimations have been studied in the visible, 
thermal infrared, and microwave parts of the electromagnetic spectrum. 

(i) Thet Photographic Method ~ The proportion of incident visible radiation that 
is reflected from soil particles decreases with increasing thickness of water film 

around the particles. This is a result of several phenomena (Planet, 1970), the most 
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Impori-ant one being multiple reflection at the oir/water interface. Successful 
utilization of this effect is hampered by the interference of soil color and soil drying 
patterns (Winkler, 1966; Cooper and Smith, 1966), negligible subsurface penetration 
of the visible radiation, atmospheric interference, etc *; consequently , only qualitative 
estimates of the surface moisture content have been made by this approach to date . 

(ii) The Infrared Method - The flux density of radiation emitted by the soil in the 
thermal infrared region is proportional to the fourth power of soil temperature. Inasmuch 
as water has an important effect on soil thermal properties, soils with various moisture 
contents may sometimes be delineated on thermal Imagery (Hunter, 1969), especially 

if sequential coverage is available. Interference of atmospheric and other variables 
also influencing the radiation recorded by an infrared sensor represents an important 
drawback with regard to obtaining quantitative moisture estimates. The depth of 
penetration is very small for thermal infrared radiation but some response to subsurface 
moisture exists, mainly due to the diffusivity variations with changing moisture 
content (Myers and Heilman, 1969; Blanchard et al,, 1974f Idso et al,, 1975), 

(iii) The Microwave Method “ The basis of the microwave method of soil water 
content determination is the difference in dielectric constant of dry soil and of water. 

While the real part of the relative dielectric constant of dry soil is less than 5, that of 
water can be higher than 80 at low microwave frequencies. Consequently, the dielectric con~ 
stant of moist soil will increase compared to that of dry soil;the increase should be in some 
proportion to the amount of water present. Figure 1 shows real and imaginary parts of the 
relative dielectric constant of loam as a ; motion of volumetric water content measured 

at three frequencies. It is apparent that the dielectric constant increases with increasing 
moisture content. However, the imaginary part responds differently to frequency 
change than the real part (Figure 1), These differences become Important in interpreting 
microwave signals in terms of soil moisture (Ulaby et al,, 1974) . 

For a homogeneous soil with a smooth surface the amount of reflected microwave 

energy is given by the power reflection coefficient p (Cihlar and Ulaby, 1974), Figure 2 

demonstrates that the reflection increases with increasing moisture in the soil. Although 

the data cover a wide ronge of frequencies and soil textures, no systematic irend of the 

power reflection coefficient with frequency or texture is apparent. If the soil is in 

addition in thermodynamic equilibrium with its environment, then emissivity e is equal 

to 1 - P ; this parameter is indicated on the second vertical axis in Figure 2, Moisture 

3 3 3 3 

content change from 0.0 cm /cm'^ to 0,45 cm /cm can be seen to cause emissivity 
(and power reflection coefficient) change of 0.4. 
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The two above parameters, p and e , ore directly related to the outputs of 
microwave sensors. For example, the scattering coefficient a° measured by radar at 
nadif over a smooth surface is given by (Peake, 1959) 

C? = 4 TT p . (1) 


Similarly, radiometric temperature T is proportional to emissivity; 

op 




( 2 ) 


Tp is physical soil temperature, and is microwave radiation emitted by the upper 
hemisphere and reflected by the surface in the direction of the radiometer. Atmospheric 
attenuation effects were neglected in Equation (1) and (2), 

Under natural conditions, soil surfaces are generally not smooth; Equation (1) must 
therefore be modified to account for surface scattering. Secondly, the soil profile is 
usually not homogeneous because of variations in moisture content, soil texture, bulk 
density, etc .; subsurface reflection and scattering may therefore occur, thus effecting 
the signal magnitude measured by the sensor. Thirdly, since the measured signal cCn con** 
sist of contributions from the surface as well as subsurface, it follows that the vertical 
distribution of soil moisture (due to reflection coefficient P and emissivity e ) and soil 
temperature (due to both radiometric temperature T and temperahire dependence of the 
dielectric constant of water ) can become important. Fourthly, the maximum depth at which 
moisture can be measured from satellite altitudes appears limited to the upper part of the 
soil profile (see Chapter 4). Consequently, if moisture information about deeper layers 
is also needed, the microwave method must be supplemented with another technique. 


1.2 PROBLEM STATEMENT 


Due to limitoHons of phofographic remofe sensing fechniques, fhe only potonHally 
useful mefhods of Iorge~scale soil moisture determinations are water balance and microwave 
and infrared remote sensing techniques. 

Studies of the microwove method of soil moisture determination conducted 
so far centered primarily on experimentally documenting the response of microwave 
sensors to soil moisture, and, more recently, on developing models which would 
allow prediction of soil moisture content from microwave signal measurements; a list of 
previous studies has been compiled by Cihlar and Ulaby (1974) . The above discussion 
suggests, however, that an understanding of problems involved in the microwave 
apporaeh to soil moisture content determination as well as the development of an 
operationally efficient algorithm for an extraction of moisture information from the 
microwave data will not be completely successful unless consideration is also given 
to moisture and temperature regimes of the soil. It is in this direction that the 

present study was undertaken, its objective being to establish a basis for incorporating 
soil moisture and temperature regimes into the microwave method of soil water content 
determination. The study was limited to bare soil'where the microwave method appears 
more accurate. In addition, studies of moisture and temperature regimes of a vegetated 
soil require an approach different from that for bore soil . Consideration of both bare 
and vegetated soil would be beyond the scope of a project of this type. Furthermore, 
the ground water level is assumed to be sufficiently deep so that it does not affect 
moisture flow within the soil profile. 

Moisture and ten^erature regimes have been the subject of research in numerous 
countries, by many investigators, and for a number of years. Rather than attempting 
to encompass all the various aspects that hove been dealt with in the past by the 
investigators involved, the approach taken here was (i) determining specific parts of 
these regimes which are of importance to the microwave method, and (ii) attempting 
to provide answers through research based on relevant data collected either for this 
study or by other investigators. The following specific problems were identifieds 
(i) Vertical moisture changes near the surface of a bare soil and 
the possibility of accurately predicting them by means of a 
mathematical model . Such a model could be used as a 
counterpart or complement of the microwave method in a 
large^area moisture estimation scheme . 
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(ii) Diurnal voriaHon of soil mois^ure near the surface. Since the 
microwave meriiod is time'^specific, large diurnal changes in 
moisture could cause error in predicting the actual water 
content for a given location. 

(iii) Vertical and temporal moisture changes at greater depths in the 

soil and the extent to which they are related to near surface moisture 
fluctuations . If such a relationship exists, subsurface moisture 
might be predicted without necessitating its direct measurements. 

(iv) "Mean maximum*' bare soil temperature variations to be expected 
at mid-latitudes as a result of varying surface slope and aspect, 
time of day, soil moisture content, and soil depth. Knowledge 
of these values may facilitate estimating errors to be expected 

in moisture determination by the microwave method due to soil 
tenr^erature, as well as ways of minimizing th^e errors. 

(v) Brightness temperature and power reflection coefficient changes 

of a bare soil in which both temperature and moisture vary diurnal ly. 

This step combines some of the above results in terms of the microwave 
method. 

(vi) An algorithm for soil water content monitoring over large areas. 

The specific problems outlined above also suggest convenient subdivisions for 
presenting the results. Chapter 2 deals with soil moisture regime (problems (i) through 
(ill)). Chapter 3 describes results of soil temperature simulations, while the coniblned 
effect of moisture and temperature as it is reflected in the brightness temperature and 
the large area moisture estimation scheme are discussed in Chapter 4. A limited 
discussion of power reflection coefficient changes is also included In Chapter 4. Chapter 
5 presents o review of experimental investigations using active and passive microwave 
sensors to determine soil moisture content. 


CHAPTER 2. 







SOIL MOISTURE REGIME 

2.1 ALTERNATIVE APPROACHES TO STUDYING SOIL MOISTURE REGIME 

Various aspects of soil moisture regime hove been studied in severol disciplines 
for a number of years. Most investigations of bare soils were made by soil physicists, 
although meteorologists and hydrologists also contributed to an understanding of soil 
water behavior, particularly under field conditions. Studies of soil moisture may be 
clossified as atomistic vs. holistic, field vs. laboratory, and empirical vs. theoretical, 
aiTiong others. This discussion is approached with the assumption that in a large”area 
determination of soil moisture, only the holistic, field, and empirical or semi*empirical 
concepts are of interest at this stage. The preference for a holistic approach is determined 
by the iact that soil water is a part of a dynamic system and is therefore affected by many, 
often interacting, parameters. The other choices (field and empirical) will be briefly 
commented upon before o review of pertinent llteroture is given. 

Laboratory experiments, although valuable for explaining the mechanisms involved 
in soil water behavior, are usually too simplified compared to the field situation. As” 
sumptions commonly mode include isotheraiol conditions, uniform profile, nonrodiotive 
regime, constant evaporation rate, etc. It should be noted, however, that some of 
these studies gave rise to simplified theoretical models of sol! water which will be 
discussed later. 

Theoretical equations can, in principle, accurately describe soil water status. 
Because of the d)oiamism of soil water, it has been argued (Nielsen et a!., 1967) that only 
numerical, as opposed to analytical, methods are of potential practical usefulness. Either 
technique involves parameters which must be evaluated for a given soil before the 
theoretical models can be used? these parameters may vory with distance, time, and soil 
condition (Nielsen et al., 1973). In addition, solution of the differential equation of 
v/oter flow requires knowledge of the boundary condition which, especially for a wet or 
moist soil, is a function of otmospheric conditions. Due to these characteristics, 
theoretical models appeor to have little usefulness for predicting field moisture changes 


10 




ah hhe presenh. Other models, developed after simplifying assumptions were made, 
have been shown to yield soil moisture estimates which agreed with experimental 
measurements (Black etal., 1969). Some authors (Gardner, 1973; Gardner and 
Gardner, 1969) have demonstrated that the simplified solutions can be used to predict 
total evaporation for a homogeneous soil. However, the practicality of these methods 
under field conditions is yet to be established. 

Apart from direct measurements, almost all traditional methods of soil water content 
determination under fie'd conditions are based on a budgeting approach. The basic 
relationship Involved in determining change in storage S is 

-§- = I (t) - O (t) (3) 

where I(t)and 0(t) refer to all inputs and outputs in the soil at time t. Equation (3) may 
be applied to a system of any size provided that l(t)and 0(t) can be determined. For a 
three dimensional body of bore soil the boundaries of which extend above the surface. 
Equation (3) can be written in the following form for a period At; 


A S P + I + F. “ F - AE, 
1 o 


(4) 


where 


S = change in soil water storage; 

P = precipitation and dew during At; 
I - irrigation during At; 

AE - actual evaporation during At. 


Assuming subsurface lateral homogeneity, the inflow F. and outflow F^ hove surface 
horizontal (F.j^, F^j^) and subsurface vertical (F.^, F^^) components. Usually runoff 
R and drainage D are defined respectively as 


and 


« = foh-fih' 


If one assumes horizontal uniformity (R = 0), Equation (4) can be written as 


(5d) 

(5b) 


AS=P + I-AE- D. 


( 6 ) 
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The amounf of waj'er held in fhe soil of fhe Mme (f +AI') can fhen be computed os 


S {t +At) = S (0 + AS 

= S (t) + P + I " AE “ D • 


(7) 


Individual components of Equation {7} can be determined in various ways , 

Initial Water Storage 

Initial water storage S(t) is determined by direct sampling or assumed to hove a 
certain value, for example after snowmelt. If the computation is performed for a 
sufficiently long time prior to the time interval of interest, the choice of S(t) becomes 
insignificant. 

Precipitation 

Precipitation amounts are measured by rain gauges and extrapolated into areas 
between gauges to obtain '*areal“ rainfall. However, there is considerable uncertainty 
associated with this extrapolation. A measure of this uncertainty can be obtained from 
dense raingauge networla. Figure 3 shows the ratio of maximum average areal rainfall in 
an area to the average maximum point rainfall over the area plotted vs. the area size. Data 
for 30 minutes, 1 hour, and 24 hours duration were measured at networks located throughout 
the United States, while curves for 3 and 6 hours were interpolated. It is apparent thot 
GS the area size increases and duration of a storm rainfall decrecQes, the differences 
between point and areal rainfall become more pronounced. For example, in a ITOsq. mi. 
area, a 30 minute rainfall may be 1 ,67 times higher at some point than the average rain- 
fall over the area (Figure 3). The location of maximum rainfall is believed to be rondom (U.S, 
Department of Commerce, 1958). Hershfield (1969) concluded from dense networks data 
collected in various parts of the U.S. that in general, the standard deviation of point 
rainfall increases with an increasing rainfall amount; that is, the coefficient of vdriobil ify 
remains relatively constant os the area size increases. Extreme variations of storm totals, 
actually observed at 20 dense networks and during storms 30 minutes to 39 hours in 
duration, are shown in Figure 4 by an enveloping line os a function of the shortest 


Figure 3i Vdridbility of raih^^ as a funotion of area and rainfqlf duration. From U ,S. Department of Commerce (1958) 








Stonn rainfall di^ (inches) 



Dlst:ance between Talngages (miles) 


Figured. Storm r^infoU differences as a function of dhtdnce rqingduges* From Hershfi eld (1969)n 






distances between raingauges. Since these data were based on a relatively small 
sample, a larger search using denser networks would provide measurements to raise 
the enveloping line (Mershfield, 1969); nevertheless. Figure 4 suggests that the rainfall 
could differ by 2.54 cm over a distance of 0.16 km* 

From the above data it is evident that precipitation can be a source of considerable 
uncertainty in using the water balance approach for soil moisture content estimation. 
Various methods have been developed for estimating areal rainfall from limited raingat^e 
measurements (e«g., Mutchinson and Walley, 1972; Chidley and Keys, 1970), but the 
predictions given by these models are limited by the input data. 

Irrigation 

In contrast to precipitation, irrigation amounts are usually well known. The 
areal distribution ts olso more uniform although, depending on the method of application, 
differences in the amount of water added to various parts of o field can occur. 

Evaporation 

Three conditions must be satisfied before evaporation from bare soil can occur. 
First, water must be available at or sufficiently near the soil surface. Secondly, latent 
heat must be available at the same sites. Thirdly, a positive water vapor pressure 
gradient must exist between the soil surface and the oir above It. 

The evaporation from an initially wet bare soil occurs in three stages (Penman, 
1941; Philip, 1957; Philip, 1967). The first stage is characterized by a high evap- 
oration rote with the maximum evaporative water loss limited by atmospheric condit .. 

The second ste^e is initiated when soil water cannot be transmitted fast enough to t ^ 
evaporation sites . Here evaporation rate decreases, surface dries rapidly, and wat 
transfer in vapor form becomes important. The third stage is typified by a low, 
relatively constant evaporation rate controlled by forces of adsorption at the liquid- 
solid interfaces (Lemon, 1956), Idso et al. (1974) demonstrated the existence of these 
three stages under field conditions. They showed that the third stage is Initiated at 
the surface water content that conrespcxids to a retention of two molecular layers around 
surface soil particles. The trcnsitions between stages 1,2 and 2, 3 could be identified 
on the basis of albedo measurements. 
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The above discussion indicai'es i'hai' bare soil evaporation is either primarily 
atmosphere*! imited or primarily soiI*limited. Accordingly, the process has been 
shidied extensively by both meteorologists and soil scientists, omong others. Numerous 
approaches with various degrees of empiricism and with more emphasis on some portion 
of the evaporative process than on others have been developed. Nunn et al . (1970) 
summarized four mass transport formulas, twelve aerodynamic formulas, an eddy correlation 
method, an energy balance method, and five empirical methods for computing evaporation. 
Equations of Baler and Robertson (1965), Christiansen (1968), and Eogleman (1967) should 
be added to the last group. These formulas require various inputs, but in general, the 
more complex ones are more accurate and can be applied to shorter time periods. The 
simplest formulas have been developed for a period of one month, although their use 
for one day intervals can give satisfactory results (McGuinness and Bordne, 1972; 

Eogleman, IS^I). 

The soil**limited evaporation stages have been dealt with at various levels ranging 
from theoretical to empirical . Many of the theoretical studies are of limited interest here 
because they do not yield a practical method for evaporation determination In the field. 

Two recent studies warrant mentioning, however. Black et al . (1969) predicted cumu* 
iative evaporation from a Plainfield sand to within 0.3 cm of water for a period of 3 months. 
The prediction was based on an analytical solution of the flow equation. Gardner (1974) 
used fractional water loss vs. square root of time curve as the basis for predicting 
evaporation from a fallow field soil . !t should be noted that in both studies, the necessary 
relationships were derived experimentally for the same soil, i.e., the "training" and 
"testing" sets overlapped; thus the predictive capability of these techniques for other 
soils has not been fully demonstrated. The studies did show, however, that simplified 
theoretical relations can be used successfully under field conditions. 

On the more empirical side, the most common approach to studying the soil** 
limited evaporation stages have been attempts to determine the relationship between the 
actual (AE) and the potential (PE) evaporation. Various relationships have been derived 
from studies of evaporation or evapotranspiration (Baler, 1967; Selirio, 1969). Selirio 
(1969) argued that the type D*reIationshlp (exponential decrease of AE when soil 
moisture Is reduced to a certain level), presented by Lemon (1956), Philip (1957), Holmes 
and Rdsertson (1960) and others, is valid for both bare and vegetated surfaces. Selirio 
and Brown (1971) obtained close correspondence between bare soil moisture measured and 
calculated using the type D*«jrve. On the basis of data from five different lysimeter 
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sl-udies, PriesHey and Taylor (1972) plowed Jhe AE/PE raHo againsi* fhe accumula^ed 
actual evaporaHon minus precipitalion difference. They concluded fhaf fhe AE/PE 
ratio equals unity until the accumulated difference reaches some critical value and then 
decreases linearly; the critical value was 0,2 cm in one bare soil case. According to 
Baier (1969), the AE/PE relationship depends on the moisture characteristics of the soil. 

Drainage 

Soil water in excess of that which can be retained for extended periods of time 
drains out of the soil profile. Its amount can be determined directly by measuring moisture 
content changes in the profile v/hen evaporation is prevented (Wilcox, 1959; Ogata and 
Richards, 1957; Hillel, 1971; etc.). Alternatively, the drainage amounts can be estimated 
from water flow theory using experimentally measured hydraulic gradient (Stone et al., 
1973), However, without direct measurements or some knowledge of soil profile charac- 
teristics, drainage estimates become difficult. Consequently, when computing moisture 
storage in the subsoil, it is often assumed that water in excess of a storage capacity of 

the soil profile drains into deeper layers. 

Various water balance studies have been made, for both bare and vegetated 
conditions . The studies involve Equation (7) and differ primarily in the way in which 
individuol components ore determined. Black et al . (1969) based the estimates on soil 
water flow theory and on empirical determination of the parameters needed. Quashu 
(1969) computed moisture depletion from soil initially at field capacity by an exponential 
decay relotion in which a "depletion coefficient” was an experimentally evaluated 
lumped parameter; this coefficient expressed a combined effect of climatic^ biotic, 
and soil variables. Eagleman (1971) calculated AE from a formula involving powers of 
PE and of moisture ratio (ratio of available water to the maximum available water in 
the soil profile). Baier and Robertson (1966), Selirio and Brown (1971), and Vanderlip and 
Brown (1974) used layered models, in which water balance is computed separately for 
a number of layers within the soil profile. In general, these studies demonstrated 
that water content within the soil profile can be monitored using the water balonce 
method, but the accuracy achieved varied with the complexity of the model and the 
input parameters used. For application to larger areas, rainfall variability presents a 
major problem. In addition, runoff enters as an important factor on sloping surfaces. 

It should be noted that hydrologic models of watersheds (Holton et al., 1974; 

■s* 

Smith and Lumb, 1967) are also water balance models. Some of these models use 
Antecedent Precipitation Index os the measure of soil moisture deficiency; Its value is 
determined as 
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API (0= API (^- 1) K + P (t). 


( 8 ) 


where API(i-) = value of API for day f, 

K = recession factor * 

As stated in section 1 .2, three specific soil moisture problems were pursued 
in this study. 

(i) Development and testing of a layered water balance model for 
predicting bore s6il moisture chonges in the upper portion of the 
soil profile. Layered models are preferrable to models which regard 
the entire soil profile os one layer because they allow prediction 
of the temporal moisture changes as a function of depth; moisture 
varies with depth in both bare and vegetated soil . However, existing 
layered water balance models do not have sufficient depth resolution. 

For example, the bare soil model of Baier et a| . (1972) consisted of 
three layers, and that of Selirioond Brown (1971) of four layers. In 
addition, layered models require depth-dependent coefficients as an 
input, and these coefficients are not readily available. 

(ii) Diurnal variation of soil moisture near the soil surface. Since the 
time increment for which moisture changes are predicted by water 
balance models is one day or more, diurnal water content variations would 
not be accounted for by this technique. Therefore, from the view- 
point of real-time measurements, the magnitude of the diurnal soil 
moisture changes under various conditions becomes important. 

(iii) Vertical and temporal moisture changes at greater depths ("subsurface") 
and the extent to which they are related to near"surface moisture 
fluctuations. This problem refers to the possibility of estimating sub- 
surface moisture without necessitating its direct measurement. 


2.2 near-surface MOISTURE REGIME 


2.2.1 Layered Water Balance Model 

2.2. 1.1 Development 

Development of a layered water balance model involves (i) specifying 
methods for computing inputs (precipitation, drainage) and losses (evaporation, 
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drainage) into individual layers, and (ii) combining these into a model. Since 
evaporation 1s the main depletion mechanism near the surface, its determination Is 
crucial for constructing an accurate model; therefore, it will be dealt with first, 
2.2,1 ,1,1 Evaporation ~ Several methods for computing evaporation have been 
described in section 2.1, Because of its simplicity, the approach based on AE/PE 
ratio as a function of soil moisture was adopted here. If the soil profile consists of 
n loyers, the amount of water lost by evaporation from a soil layer ] on day i, AE.., 
can be calculated as (Baier et dl., 1972) 


AE.. = k. C.. PE. , 
«J J 'J ' 


( 9 ) 


,th 


where C.. = coefficient accounting for those soil characteristics in the j layer 
which affect evaporation loss; 

PEj = potential evaporation on day i; 
k. = depth coefficient for the loyer. 

The total amount lost from n layers, AE., Is equal to 

n 


AE. = PE. Xk.C.. , 
' ' J 'J 


( 10 ) 


If moisture content is sufficiently high, AE. = PE. (section 2.1), and C.j = 1 
(Baier etal., 1972), Therefore, 


S'." 


(II) 


Equation (11) shows that coefficients k. distribute the total evaporated amount AE. 
to individual layers within the profile. The values of k. can thus be determined 
directly from Equation (9) when C.j = 1 , When actual evaporation falls below potential 
evaporation, C.j becomes less than one and its value will be a function of soil drying 
characteristics (Section 2.1). 


I ill, 
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In order to use Equation (9) , one must obtain estimates of k. and of C.. . 

As noted before, k. values published in the literature (Baier et al., 1972; Selirio 
and Brown, 1971; Vanderlip and Brown, 1974) do not provide the spatial refine" 
ment needed . It therefore became necessary to obtain new kj estimates with 
adequate depth resolution. Since this method of AE| computation is essentially 
empirical , the new estimates must be derived from accurate measurements in order 
for Equation (9) to yield adequate results . 

In the course of an experiment dealing with near surface water fluxes in a 
bare soil, researchers at the U.S .Water Conservation Laboratory in Phoenix, Ariz(xia, 
acquired a set of detailed measurements of soil moisture. These included volumetric 
soil water contents m. for seven depth increments (0-1 cm, 1-2 cm, 2-3 cm, 3-4 cm, 

4-5 cm, 5-7 cm, and 7-9 cm) at 20 minute intervals for a nun^er of days, and total 
actual evaporation AE. (see Appendix B for their experimental procedure). Jackson 
(1973) presented daily values of moisture for a 37-day experiment conducted in March 
and April of 1971; the some data were provided for the present study in tabular 
form*. Since these values were computed from individual measurements between 
0000 (midnight) and 2400 for any given day, the averages were considered representative 
of conditions at 1200* hours . Values of were computed as the difference between 
m . for adjacent days, and corresponding m. was determined as the average of m.'s 
for adjacent days. Pan evaporation measurements were obtained from the University 
of Arizona Citrus Experimental Station (latitude 33°23'N, longitude 111^58’W for the 
same time period. Under the conditions of the evaporation experiment (Appendix B), 
the potential evaporation (PE) could be expected to occur at the levels close to pan 
evaporation because of the large thermal advection. This was also confirmed by the 
ratios of actual and pan evaporation which were close to 1 .0 for the first three days 
in the March, 1971, experiment. Consequently, the measured pan evaporation values 
were considered equal to potential evaporation (PE) for all data collected at the U. S. 
Water Conservation Loboratory. Since the PE values were recorded at 0800 hrs,, the 
measured PE values were partitioned for the 1200 to 1200 periods. Consequently, 
sets of m.., AEj, were available for an evaporation period of 37 

days. 


*Courtesy of Dr. R, D. Jackson, U. S, Water Conservation Laboratory, Phoenix, Ariz, 
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To determine values of k., ratios AE.j/PE. were plotted against the corresponding 
'”i/* for eoch of the seven layers between 0 and 9 cm. Data points from a July, 1970, 
experiment were plotted in Figure 5 and others for comparison only; it should be stressed 
that the March, 1971, data alone were used in deriving k. and X coefficients. Figure 5 
shows the AE.j/PE. ratios for depth 0“1 cm: at the highest moisture content, the value is 
approximately 0.13. In this fashion, preliminary estimates of k. were obtained for 
depths 0 to 9 cm. These estimates were preliminary because even for the highest moisture 
content for which data were available, AE. was less than PE., their rotio being 0,90. 

For this reason, a second k. set was obtained after dividing the preliminary estimates by 
0.90. 


In the computations described so far, no account was taken of the layers below 
9 cm from which evaporation also occurs. The total thickness of the evaporation zone 
in bare soil is not constant but depends on soil water diffusivity, moisture gradients, 
evaporation demand, etc. Data from the same experiment indicated that the plane 
of zero water flux varied between 15 and 35 cm (Jackson et al ., 1973) . Laboratory and 
field experimental data hove been presented which show that the bulk of evaporated 
water originates in the 0“15 cm to 0"70cm layer (Benoit and Kirkham, 1963; Keen, 1922; 
Fritton et al,, \ 970 } Hanks and Gardner, 1965; Willis and Bond, 1971 ), In this study, 
a constant thickness of 30 cm was chosen; this also corresponds to an upper moisture 
zone thickness in some hydrological models (e.g.. Smith and Lumb, 1967). Since no 
detailed measurements were available for the layers between 9 and 30 cm depth, it was 
assumed that k. values decrease linearly in that region. These extrapolated values, 
together with the experimentally derived values, are shown in Figure 6 and Table 1 , 

Given k.*s from Table 1, values of C.. (Equation (9)) can be computed for all 
J J 

measured moisture contents; for any day i and layer j. 


C.. 

•J 


ae.. 

LL_ 

k.PE. 


(12) 


C.j will thus decrease from unity when AE. = PE. to zero when AE. = 0. Figure 7 and 
8 show C|. values calculated for data from several experiments (all on Adelanto loam) 
at two depths, 0-1 cm and 2-3 cm. The volues of C.. remained constont to some value of 
moisture content and then sharply decreased; this suggests that the exponential decay 
relation between C.. and m. (tjqse D"curve, section 2.1) holds for these data. There 
are two noteworthy regularities, however. First, the critical moisture content below 
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Ratio of Actual (0-1 c 





Figure 5. Relationship between the ratio of actual (depth 0~ 1 cm) and 
potential evaporation and soil water content. 


lAr- 
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Depth j, cm 



Figure 6. Coefficient k as a function of depth. 
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Coefficient C;;, Dimensionless 






which C.j decreased to less than one, is lower near the surface than at greater depths. 
Secondly, the decrease of C.. becomes more rapid os the depth increases (Figure 7 and 
8 ). ^ 

The exponential decay relationship may be expressed as 


C.. = 1 if m.. 10.95 SC. 03a) 

U 'J J 

m.. ” 0.95SC . 

= X otherwise. (13b) 


where SC. = storage capacity of the layer. Selirio (1969) used X = 127.42. 
Equation (13b) represents a family of curves, however, which differ primarily 
in the rates of decoy. Attempts to colculate X from the experimental data were not 
success^ I because of the large scatter of dota points (Figure 7 and 8). Therefore, 
values of X were derived by trial and error so as to represent the scatter of C.. points 
for each layer. It may be noted that Equation (I3b) can also be written os 


m.. “ 0.95SC. 

C.. = exp ( a — ^ ) , (14) 

a = InX. 

Values of X and a for 30 layers 1-cm thick are given in Table 2. Values 
for depths greater than 9 cm were extrapolated on the basis that os X increases by a 
factor of 2, so does the thickness of the soil layer to which this X value applies. 

2.2.1 .1 .2 Distribution of Precipitation ~ The technique used to distribute precipitation 
to various layers was identical to that of Baier et ol (1972). As water infiltrates below 
the surface, each layer is saturated to its storage capacity before water drains through 
it into the next layer; water remaining after saturating the deepest layer was considered 
os drained out of the profile and did not enter any computations. 
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2.2,1 .1.3 Model ~ A flowchart' for the wafer balance model is shown in Figure 9, 

The actual compufafion of soil moisture change for a given day i consisted of three steps: 


0 ) 


(ii) 


(Hi) 


Assignment of depth coefficients kj to each layer j . For a 

sufficiently moist soil, kj's were equal to those in Table 1 . 

If the near-surface layers dried out below a "critical water 

content" (CRWC), k. values were shifted down one layer for 

each "dry" layer. This modification was based on the 

observation that after moisture content of the 0-1 cm layer 

3 3 

decreased below approximately 0,06 to 0,08 cm /cm , 

moisture loss from the 1-2 cm layer increased. According to Jackson 

3 3 

(1973), moisture content of 0.06 cm /cm for Adelanto loam 

corresponds to approximately two molecular layers, and that of 0.08 
3 3 

cm /cm is equivalent to a tension of 200 bars. In the following 

3 3 

sections, CRWC was assumed equal to 0.08 cm /cm for loam 

3 3 

and silt loam, and 0.10 cm /cm for clay loam. 

Computation of actual evaporation AE.. and total evaporation 

AE. for NL layers. This confutation followed Equations (9) 

and (13). Values of AE.. were set equal to 0.0 if m.. wos less than 
q ^ q U 'J 

0,03 cmycm'^. 

Distribution of precipitation measured for day i was carried out as 
previously explained (section 2. 2 *1.1. 2; Figure 9) with one 
modification: the deepest layers reached by the rainfall were not 
recharged completely so that a more realistic moisture gradient 
with depth could be obtained . 


2,2,1 .2 Testing 

The water balance model described in previous sections wos tested on four sets 
of data (Table 3). The first set (TSl) was used to derive values of k. and C... The 
second set (TS2) was obtained in an experiment carried out in July, 1970 , which was 
otherwise similar to the March, 1971, experiment (Appendix B), The third (TS3) and 
the fomth (TS4) testing sets were obtained during measurements of radar backscatter 
from bare soils near Eudora, Kansas; data for the 1973 season (TS3) were presented by 
Cihlar (1973), For these comparisons (Tables 5 though 8) predicted values were averaged 
so as to coincide with the depth increments sompled in various testing sets. The follow- 


29 




WGINAL PAGE 18 
OP POOR QUAU3Y 






Table 3. Paramel’ers of four dai’O sets (TS1 through TS4) used to test the layered 
water balance model . 


Data 


Date 

Number of Days 


Potential Evaporation 



Location 

and 








1 



Set* 




When 

Soil 

Samples 

Taken 

Depth Increments 


Statistics (cm) 



No. 

Soil 

From 

To 

Totol 

sampled (cm) 

Source 

Mean 

Standard 

Deviation 




1 

Phoenix, Az. 

W71 

4/Q/71 

33 

15 

0-1, 1-2, 2-3, 3-4, 4^5, 

Pan 

0,526 

0.069 




Adelanto loam 





5-7, 7-9 

evapora- 

•k 

Hon 




2 

Phoenix, Az . 

7/12/70 

7/17/70 

6 

6 

as above 

Pan 

1.077 

0.063 

1 



Adelanto loom 






evoporo" 

, . * 

tion 




3 

Eudoro, Ks. 

7 / 25/73 

9/5/73 

43 

8 

0-1, 1-2, 2-5, 5-9, 9-15, 

* 

Computed 

0.560 

0.104 



Eudora silt 
loom 





15-25, 25-35 





4 

Eudoro, Ks. 

Pawnee clay 
loom 

e/25/72 

9/15/72 

28 

3 

0-2, 2-5, 5-9, 9-15 

Computed 

0.445 

0.095 



5; 




*See text for explanations 





















ing data sources were employed to calculate potential evoporation for TS3 and TS4. 
Temperatures, relative humidities, wind speed, and rainfall were either measured at 
the experimental site or computed from linear regression equations developed on the 
basis of on~site measurements (dependent variable) ond of the University of Kansas 
meteorological station measurements (independent variable) for the overlapping 
periods. Estimates of percent sunshine were obtained by overagir^ values reported by 
the U. S. Weother Bureau Airport Stations in Topeka and Kansas City, respectively. 
Values of daylight duration and fotal radiation outside the atmosphere were taken 
from List (1951). Using these data, potential evaporation values were computed by 
methods of Penman (1948), Baler and Robertson (1965), Christiansen (1968), Eagleman 
(1967), and Thornthwaite (1948); however, only results obtained from the Penman 
equation were used in the water balance calculations. 

Figure 10 shows computed and measured (TSl) moisture contents for three layers 
(0-1 cm, 3^4 cm, 7“9 cm). The predicted values were in good agreement with meosured 
values, especially for the period 4 to 16 days after irrigation. For the remaining two 
measurements (23 and 37 da)re after irrigation), soil at greater depths dried out faster 
than the model predicted. This suggests that the "sinking" of k. values for lower moisture 
contents did not raise AE.. values sufficiently. It is possible that the CRWC value should 
also increase with depth, thus allowing for faster "sinking" of k. coefficients. The 
difference between computed and measured values was worst at the depth of 3“4 cm 
as indicated in Table 4. These data show that 23 days after irrigation, prediction was 
satisfactory for all depths except the 3“4 cm layer, while its accuracy decreased for 
the adjacent layers during the subsequent period. 

Table 5 gives a statistical description of the correspondence between measured 
and computed values for TSl; all pairs of values ore plotted in Figure 11 . The good 
predictive copabiUty of the model is indicated by the high correlation coefficient 
( S 0.959 for individual layers) and low stmdard error of estimote ( ^ 0,0111 for 
individual layers). An interesting ond repeatedly occurring feature of the predicted 
moisture contents was that the slope of regression lines for individual depths was less 
than unity; that is, the model overestimated actual evaporation for a moist soil but 
underestimated it when the soil became dry. This Is probably again related to k. 
values although the slope was also below unity for the surface layer in two of the 
four testing sets (TSl, TS3; see Tobie 5,7). 
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Table 4. Differences between com^Hifed and measured soil water contents 
at 16, 23, and 37 days after irrigation, data set TSl . 


Day 

Since 

Irrigation 


Difference between computed and measured moisture content for 
depth increment (cmVcm^). 


- 0.010 

- 0.011 


+ 0.010 

+ 0.013 


0.000 -0.015 - 0.009 


+ 0.011 

+0.038 

+0.036 


Table 5. Statistics for the layered water balance model performance 
data set TSl. 


Depth 


Mean 

Measured 

Computed 

0.094 

0.092 

0.148 

0.142 

0.183 

0.184 

0.199 


0.208 


0.217 

0.223 

0.224 

0.228 


Slope 


.936 
.887 
0.914 
0.717 
0.774 


Correia" Standard 
Intercept tion Coef- Error of 
ficient JEstimate 


0.004 

0.010 

0.017 

0.065 

0.052 

0.054 

0.028 


.984 
.986 
0.975 
0.959 
0.959 
0.965 
0.976 


0 . 

0 . 

0 . 
0.0073 
0.0080 
0.0073 
0.0061 


0.182 0.186 0.940 0.015 0.976 0,0135 




























Measured Soil Water Content, cm^/cm^ 


Figure 11. Comparison csiF computed and measured soil water contents. Dote 
Set TS1 . 
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The correspondence between moisture values for Adelanto loam collected In 

July, 1970 (Jackson, 1973) and predicted values is shown in Figure 12. Except during 

two days for the 1-2 cm layer, the pairs of values agreed closely. Since the measurements 

were made only for the initial evaporation period, the accuracy of prediction was 

excellent in terms of both correlation and standard error (Table 6). 

Both TS1 ond TS2 were related to the model derivation, so that a close 

correspondence between computed and measured moistures may not be surprising . 

Indeed, a cursory exominotion of the statistics for individual depths of TS3 (Table 7) 

reveals substantial decrease in accuracy . However, measured moisture values used for 

the statistical comparison with computed values were means of one or more samples, 

and the individual samples exhibited considerable variations around the means. The 

sampled mean values were plotted against the computed ones in Figure 13. In oddition, 

error bars were placed on the extreme measured values when more than one sample were 

taken . If cases with only one sample are ignored and assuming for the remaining points 

that the actual soil moisture could be anywhere within the error bars, then the maximum 

3 3 

difference between computed and measured values is less than 0.035 cm /cm for all 

Q -a 3 3 

(50) points, less than 0.025 cm/crn for 92% of all points, and less than 0.020 cm /cm 

for 84% of all points. In other words, although the correspondence between measured 

and computed values was worse than for TS1 and TS2, it was nevertheless rather close. 

This conclusion is further supported by comparing the computed values for depths of 

10 cm and 25 cm with values meosured tensiometrically (these data are discussed in 

section 2.4, and the experimental procedure used in their acquisition is described in 

Appendix B). For the period between 8/5/73 and 11/3/73, the absolute difference 

3 3 

between computed and measured values was generally less than 0,020 cm /cm . 

Con^arison between computed and measured moisture contents for TS4 yielded 

similar results: the intial scatter of points around the 1: 1 line was reduced to less than 
3 3 

0,01 cm /cm when the error bars were taken into consideration (Figure 14). Each data 
point in TS4 was a mean of two samples. The occasionally large difference between 
their moisture contents was attributed by Ulaby (1974) to the presence of a ditch which 
facilitated drainage of one site only. In addition, heterogeneity of the soil infiltration 
rates could be involved as the differences were greatest following a large rainfall. 

A statistical analysis of the predictive accuracy for TS4 (Table 8) shows good agree- 
ment between measured and cotrouted values. It should be noted that for both TS3 and 
TS4, the statistics (Table 7,8) are based on the mean measured values only; obviously, 
the parameter.' would change somewhat if the error bars were taken into consideration. 
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Table 6. StaHsHcs for fhe layered wol'er balance model 
performance, dara sef TS2. 


Depth 

Mean 

Slope 

Intercept 

Correla- 
tion Coef‘ 
ficient 

Stondard 
Error of 
Estimate 

Measured 

Computed 

0-1 

0.071 

0.065 

1.131 

-0.016 

0.974 

0.0093 

1-2 

0.141 

0.125 

0.635 

0.036 

0,994 

0.0035 

2-3 

0.183 

0.179 

0.839 

0,026 

0.999 

0.0012 

3-4 

0.201 

0.204 

0.755 

0.052 

1.000 

0.0007 

4-5 

0.208 

0.210 

0.853 

0.033 

0.998 

0.0016 

5-7 

0.214 

0.220 

0.859 

0.037 

0.995 

0.0021 

7-9 

0.219 

0.226 

0.897 

0.030 

0.997 

0.0015 

All Data 

0.177 

0.176 

1.043 

-0.009 

0.985 

0.0105 


Table 7. StaHsMcs for t'he layered water balance mode! 
performance, data set TS3. 


Depth 

Mean 

Slope 

Intercept 

Correla- 
tion Coef" 
ficient 

Standard 
Error of 
Estimate 

Measured 

Computed 

0-1 

0.187 

0.174 

0.941 

-0.002 

0.975 

0.0213 

1-2 

0.201 

0.197 

1.080 

-0.020 

0.905 

0.0319 

2-5 

0.222 

0.245 

1.441 

-0.076 

0.947 

0.0199 

5-9 

0.234 

0.252 

0.881 

0.046 

0.664 

0.0303 

9-15 

0.264 

0.257 

0.670 

0.080 

0.467 

0.0322 

15-25 

0.301 

0.283 

0.334 

0.182 

0.286 

0.0243 


0.312 

0.293 

0.463 

0.149 

0.533 

0.0090 

All Data 

0.246 

0.243 

0.913 

0.018 

0.900 

0.0282 
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Eudora Silt Loam 
1973 
Depth: 
o 0-1 cm 
A l-2cm 
o 2-5 cm 

• 5 -9 cm 
▼ 9-15 cm 
■ 15 -25 cm 

* 25-35 cm 


ulii 


ai a2 

Measured Soil Water Content, cm^/cm^ 


Figure 13. Comparison of conr^uted and measured soil water contents. Doto set TS3. 

Horizontal bars designate standard deviation, vertical bars indicote that 
only one sample was tcdcen. 
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Computed Soil Wafer Content, cm^/cm^ 


0.45 


0.0 V 02 a? 0.4 0.457 

Measured Soil Water Content, cm^/cm^ 


Figure 14. Comparison of ccHopufed and measured soil water contents. Data set 
TS4. Horizontal bars designate standard deviation . 
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The above result indicate iha( soil moisture in the 0 to 30 cm zone may be 

3 3 

predicted for 1-cm depth increments within 0,02 cm /cm . Precipitation, daily 
potential evaporation, soil storc^e capacity, and initial soil moisture contents are 
required os input. Two remarks concerning the storage capacities are in order. First, 
the appropriate value of storage capacity appears to correspond to a tension of opproxi* 
mately 0.09 bars. The moisture values used for TSl, TS2 and TS3 (Table 9) represent 
o tension of 0,09 bars, according to moisture characteristics given by Jackson (1973) 
and those obtained by measurements (Appendix B). The storage capacity initially used 
for Pawnee clay loam (TS4) was the field capacity value for a clay loam given by 
Salter and Williams (1965), multiplied by empirically determined bulk density (Ulaby 
et al,, 1974). The predicted values were too low, however; the results given In Table 
8 were obtained only ofter the above storage capacities were multiplied by 1.10. The 
low tension values cannot be explained by the necessity to provide storage space for 
precipitated water during the fiist two days after rain because no water was added in 
TSl and TS2. Although these results are consistent, further testing with data from widely 
different conditions should be made before a final value of storage capacity suitable 
for this model is accepted. It is also conceivable that the storage capacity varies 
within a moisture tension rarige similarly as the field capacity. Secondly, the storage 
capacity for the top 2 cm was 0,9 that of what it would be if the layers were located 
deeper below the surface. This adjustment, derived from Jackson's (1973) data , was 
due to the observation that for the top two layers, AE.. values were high even several 
days after Irrigation when their moisture content was much below the field capacity (Figure 
5). The close correspondence between computed and measured moisture values for these 
layers (Table 5 through 8) suggest that the trend of higher evaporation rates from near" 
surfoce layers was real; it may be ckie to the evaporative demand of the atmosphere which 


is most effective at the surface. 

Although the water balance model was constructed for 24-hour time increments, 
it was of interest to determine the accuracy with which hourly soil moisture changes 
could be predicted. If the model performs well for hourly time intervals, it could be 
used to simulate diurnol changes of soil moisture. Results described In this paragraph 
hove actually been computed using the soil temperature simuiatitxi model (Chapter 3), 
which included the water balance model for computing hourly profile moisture changes; 
they are presented here to provide a more comprehensive picture of the model's 
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Tobie 9. Storage capacities for the four testing data sets 


Depth 

(cm) 

Storage Capacity (cm^/cm^) 



TS1JS2 

TS3 

TS4 

ueprh 

(cm) 

TS3 

TS4 

0-1 

0.270 

0.275 

0.320 

9-10 

0.305 

0.423 

1-2 

0.270 

0.275 

0.320 

10-11 

0.305 

0.426 

2-3 

0.315 

0.305 

0.356 

11-12 

0.306 

0.428 

3-4 

0.315 

0.305 

0.366 

12-13 

0.306 

0.428 

4-5 

0.315 

0.305 

0.378 

13-14 

0.307 

0.428 

5-6 

0.315 

0.305 

. 0.387 

14-15 

0.307 

0.428 

6-7 

0.315 

0.305 

0.396 

15-16 

0.308 

0.428 

7-8 

0.315 

0.305 

0.405 

16-17 

0.310 

0.426 

8-9 

0.315 

0.305 

0.417 

17-18 

0.312 

0.423 





18-19 

4 

0.315 

0.420 





19-20 

0.320 

0.416 





20-21 

0.323 

0.410 





21-22 i 

0.326 

0.’408 





22-23 

0.329 

0.405 





23-24 

0.332 

0.403 





24-25 

0.335 

0.403 





25-26 

0.338 

0.406 





26-27 

0.341 . 

0.408 





27-28 

0.344 ! 

0.410 





28-29 

0.347 

0.413 





29-30 

0.350 

0.417 ’ 
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performance. The computed hourly soil water contents for Adelanto loam were 
compared with values measured on 9/20/73 (high moisture) and 10/^73 (low 
moisture content); the measured data are described in Sections 3.1 and B.2. In 
the water balance model, potential evaporation was calculated as LE/L, where LE 
is given in Appendix A and L is the latent heat of voporization. Diurnal moisture 
content changes at two depths are shown in Figure 15. These results indicate that the 
layered evaporation model predicted moisture losses rather well even for small time 
periods. Since the evaporation mode! did not compute soil water movement, the 
recharge during the night could not be predicted. It should be noted that the evap- 
oration process was simplified in that water was assumed to evaporate directly from 
a given depth rather than to flow to the evaporation sites which may be confined to 
a narrow zone in the soil (Hillel, 1971). 

It should be noted that the values of k. and C.. (Table 1, Equation 13) were 
developed and the performance of the woter balance model was tested primarily for 
high potential evaporation values . The lost two columns of Table 3 show that data 
set TSl, which served os the basis for developing kj”s and Cj.'s, had approximately 
the same evaporation levels as TS3 and somewhat higher than TS4. Since k. and C.. 
coefficients have been derived empirically, their validity for other potential evapo- 
ration levels should not be assumed without appropriate testing. The excellent pre- 
diction of water contents for TS4 where potential evoporation was almost twice as 
much as for TSl suggests that the k.'s and C..'s derived here may be valid when 
potential evaporation exceeds 0.5 cn/day. When the atmospheric evaporative 
demand is low, however, less water would be lost from the soil, and consequently 
even lower rates of water flow (associated with lower moisture contents for a given 
soil) would entirely satisfy the evaporative demand. Thus, at low potential evapo- 
ration values, one would expect that the critical storage capacity fraction would 
decrease from the 0.95 used in Equation 13. Secondly, the low rates of water loss 
might lead to a slower decrease of k. values with depth and consequently to a more 
uniform profile soil moisture distribution. These hypotheses should be tested to 
determine whether modifications of the loyered water balance model for low potential 
evaporation levels are necessary. From the point of view of water loss to the 
atmosphere, however, high evaporation conditions (where the present model appears 
to perform satisfactorily) are of greater importance. 
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Time Of Day, Hours 


Figure 15. Comparison of computed and measured soil water contents for wet and dry soils 
calculated using the layered water balance model. 



2,2,2 Diurnal Soil Moisture Variations 


Changes in soil wafer confenf dealf with in the previous sections were 
overage changes that occurred over a period of 24 hours. In this case, diurnal 
fluctuations of soil temperature may be ignored in computing evaporatic«i (Gardner 
and Gardner, 1969) because of the periodicity of temperature variations. In reality, 
however, soil moisture may change considerably during the diurnal cycle, mainly 
near the surface where the temperature gradients are at maximum. 

Field studies of diurnal moisture variations have been rather limited, in port 
due to the difficulties involved tn accurately determining small moisture changes over 
short time periods. In fact, only two studies have been reported in the literature, 
namely by Rose (1968 a,b) and by Jackson (1973) and Jackson et a|. (1973, 1974), 

Data collected during the letter series of experiments form the basis for the following 
discussion; the experimentol procedure Is briefly described in Appendix 6. These dato 
have been analyzed by Jackson and his co~workers in terms of woter flux; here we are 
primarily interested in absolute moisture levels. 

Figure 16 shows hourly moisture content changes at three depths of wet soil 
over a 24 hour period. Water content remained unchanged until sunrise and 
then began to decrease rapidly. The moisture loss started earlier and proceeded ot a 
faster rate at locations closer to the surface . Jackson et al . (1973) attributed the rapid 
surface drying during morning hours to bidirectional movement of woter: soil water from 
the surface to obout 1 cm moved upward, while that below 1 cm moved downward. In 
the late afternoon, moisture recharge from deeper layers exceeded evapwative loss and 
the water content increased again, although it did not reach the morning level. Thus 
there was a net loss in all three layers shown in Figure 1<^ this net loss was approxi- 
mated by the layered evaporation model (section 2.2.1). As the soil dried out, the 
amplitude of diurnal moisture changes was dampened (Figure 17) and the net loss at any 
depth during 24 hours decreased; os before, the amplitude was at maximum near the 
surface. Furthermore, the recharge of the near surface layer continued until sunrise. 

Data in Figure 16 and 17 represent relatively simple patterns of moisture change: | 

moisture is highest before down, then decreases until late afternoon, and subsequently 
increases again until It levels off at some pre-dawn value. In general, however, these 
changes are not so regular, especially at greater depths. Figure 18 shows diurnal 

water content changes at five different doys for a drying Adelanto loam. It is apparent I 

that the above described pattern holds well for the top 1 cm layer during the entire I 
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period of meosuremenls . For the greofer depths, the main pattern of moisture changes 
is poorly defined, partly due to the presence of higher frequency perturbations super- 
imposed on the moisture vs. time curves. When the moisture content at a given depth 

3 3 

decreased below approximately 0,120 cm /cm , however, the regular pattern observed 
in Figure 16 and 17 emerged, and the perturbations became minor (Figure 18j|; this 
condition was reached prior to day 16 (depth 1-2 cm), day 25 (2-3 cm), and day 37 
(5-7 cm), respectively. 

The above results indi'’ ’’e that the total moisture exchange (summed absolute 
values of moisture gains and losses) generally exceeds the effective moisture loss, i.e, 
the difference between final (2400 hours) and initial (0000) moisture content; the 
second quantity is equivalent to actual evaporation AE.. (section 2.2.1). In Figure 19, 
pairs of total moisture exchange and effective loss points for Adelanto loam are plotted 
agolnst the average daily moisture content; the values were determined from meosurements 
collected during seven days in two experiments (March, 1971; September, 1973), Total 
exchange was calculated for four time periods: midnight to sunrise, sunrise to noon, 
noon to sunset, and sunset to midnight. Both effective loss and total exchange were 
divided by the actual evaporation AE. to remove the effect of varying evaporative 
demand of the atmosphere. At high moisture contents, total exchange was only slightly 
greater than the effective loss. As the moisture content decreased, effective loss 
decreased to zero while the total exchange increased to more than half of the actual 
evaporation AE.. When the effective loss for near surface layers is zero, the total 
exchange should equal unity. This is not the case in Figure 19 because the 24-hour period 
was subdivided into only four increments for which the total exchange was determined; 
thus the total exchange values should actually be higher than shown in Figure 19. 

The negative loss values (i.e., moisture increase over a 24-hour period) were measured 
on day 16 (3 values) and day 37 (1 value); they could be real or due to errors in reodir^ 
off moisture contents from Figure 18, as very small quantities of water were involved. 

Figures 16 through 19 show that the diurnal change in moisture can be considerable 
depending on the depth and mean moisture level at that depth. To obtain some idea 
obout the ranges of values involved, the average daily moisture was plotted against the 
maximum and the minimum moisture contents during that doy at a given depth (Figure 20)* 
These data were taken from the same sources as those for Figure 19 and thus represent 
reasonably extreme conditions in terms of o large atmospheric evaporative demand. They 
may therefore serve os an indication of an upper limit of diurnal moisture fluctuations 
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Adelanto Loam 
Date: March 1971 
lQ/2/73 
9/2(y73 

Effective Total 
Depth Loss Exchange 

0- 1 cm 

1- 2cm 
3-4cm 
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0.20 


0.00 


o 
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° Jo ri® ■ 
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0.20L- 

0.0 


0.1 0.2 0.3 

Mean Daily Soil Water Content, cm^/cm^ 


Figure 19. Diurnal effective toss and total exchange of soil water as functions 
of mean daily soil water content. 
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Maximum and Minimum Daily Soil Water Content, cm^/cm® 


Adelanto Loam 
Date: 

March 1971 
9/20/73 
lQ/2/73 
Depth; 
o 0 cm 
a 0-0.5 cm 
• 0-1 (or 1) cm 
■ 2-3 (or 2) cm 
▼ 4-5 (or 4) cm 


1 0.2 
Mean Daily Soil Water Content, cm^m^ 


Figure 20. Maximum and minimum daily soil water content versus the mean 
doily water content at several depths. 



near (he surface of a bare soil. The variaHon was laiige for depths smaller than 1 cm 

but decreased rapidly at greater depths; for the 0~0.5 cm layer, the range was os 
3 3 

high as 0.135 cm /cm'’. Daily moisture ranges for the same data are plotted in 
Figure 21 . Except For one anomalously low value (0 cm; 9/20/73), the trend was 
consistent: maximum rarige of moistures was observed at or near the surface; (ii)with 

decrecsing soil moisture, the range decreased faster at greater depths than at the surface; 

and (iii) the range became small for depths greater than 4 cm. Typical moximum values 

3 3 3 3 

of the diurnal moisture range appeared to be 0.125 ci. /cm for (h0,5 cm, 0.06 cm'ycm 

3 3 

for 0-1 cm, and 0.02 cm /cm for a depth of 4 to 5 cm. 


2.3 SUBSURFACE MOISTURE REGIME 


The layered water balance model developed and tested in section 2.2.1 accounted 
for soil moisture changes in the top 30 cm only. There were severol reasons for this 
restriction. First, it was assumed that the dynamics of the moisture regime of a bare soil 
is greatly attenuated at depths greater than 30 cm. Secondly, soil water present near 
the surface has a more direct effect on microwave radiation as well as on the exchange 
processes between atmosphere and soil, both mass (e.g., infiltratimrt, evaporation) and 
energy (latent heot). Thirdly , the depth-of-penetration limitations of microwave sensors 
render the upper part of the soil profile most suitable for microwave remote sensing. 
Nevertheless, there are activities in which moisture content of the deeper layers is of 
importance; examples are crop production and ground water recharge. For such purposes, 
a method of soil moisture determination over large areas would be expected to yield 
information about water content status of deeper layers, in addition to those for the 
upper portion of the soil profile. 

The purpose of this section is twofold: (i) to determine whether substantial 
differences in soil moisture dynamics exist between 0*30 cm layer Qiereafter "upper 
zone") and the layer below 30 cm ("lower zone"), thereby justifying their separation 
in the previous sections; (ii) to exomlne the main features of the subsurface moisture 
regime and establish a relation between moisture regimes of the two zones which could* 
serve for estimation of moisture content of the subsurface zone without direct sampling. 



Adelanto Loam 
Date: 

March 1971 
9/20/73 
10/2/73 
Depth: 
o 0 cm 
□ 0-0.5 cm 
• 0-lcm 
■ 2-3 cm 
▼ 4-5cm 


0-0.5 cm 


Mean Dally Soli Water Content, cm^/cm' 


Figure 21. Maximum daily soil water content range versus the mean daily soil 
water content at several depths. 


5 . 


The following discussion is based primarily on field soil water content measure" 

ments. The experiment during which the data were collected is described in Appendix 

B. The data represent one soil only, Eudora silt loam. While two complete sets of 

tensiometers (SI and S2) were installed in this soil, the doto were not averaged because 

it was found that small differences existed between SI and S2 which could obscure 

subtle trends otherwise observable in both sets. Figure 22 shows the relationship between 

moisture contents simultaneously determined by the gravimetric and tensiometric methodb. 

3 3 

From the total of 32 measurements^ 16 were within 0.01 cm /cm from the 1: 1 line, 

3 3 3 3 

11 within 0.02 cm /cm , and 5 within 0.03 cm /cm . Since most of the gravimetric 
roeosurements Include some uncertainty (the upper zone data hove been discussed in 
section 2.2.1 .2), the correspondence between the two methods of measurements can 
be considered satisfactory. More importantly, the tensiometricolly acquired moisture 
measurements appeor to have accurately responded to small changes in profile moisture 
conditions (refer to Figure 28). 

2.3,1 Soil Moisture Profile Dynamics 

Changes in soil moisture content occur as a result of flow, and flow in the soil 
is proportionol to hydraulic conductivity and potential gradient according to Equation 
(15) for one~dimensional flow 

^ , 05 ) 

3 3 

where 0 = moisture content, in cm /cm ; 

3 2 

K($) = hydraulic conductivity, in cm /cm /min; 
z = ®stanw^om'tfe reference level, in crtv 
t = time in minutes. 

Since the potential gradient di/dz is determined primarily by the differences in matric 
suction, changes of tensiometric values (which are a measure of metric suction) may be 
used as an indicator of the potential of soil moisture to change. Figure 23 shows the 
highest (Q/3C/73) and lowest (9/27/73) values of moisture tension recorded during the 
period of measurements. The highest values were recorded after two weeks practically 
without rain during which the potential evaporation averaged 0.566 cn/doy. The 
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Eudora Silt Loam 
Tensiometer Sets SI, S2 
1973 
Depth: 

• 10 cm 
V 25 cm 
■ Others 


'Gravimetric” Soil Water Content, cmVcm^ 


Figure 22* Relationship between soil moistures detennined by gravimetric and tensiometric 
methods. 








lowesf- values were measured after five days of precip^^a^ion with a f'ol'al rainfall of 
17.27 cm (Figure 24). For fhe dry condition, tension was 386 cm at the depth of 10 cm, 
but decreased rapidly to 40 cm and more slowly at greater depths. Following the large 
rainfalls, tensions were low between 10 cm and 40 cm, and then increased slowly. 

The reason for increase in moisture tension at 100 cm depth is not clear, especially 
because it occurred consistently throughout the period of measurements , Differences 
between the maximum and the minimum recorded tensions are plotted as the third curve in 
Figure 23A. The range of tensions encountered during the experiment decreased 
monotonically with depth, although changes in slopes of the difference curve can be 
noted at 25 cm and 85 cm. 

Moisture content change in time is a function of both potential grodient 

(and therefore moisture tension differences) and hydraulic conductivity. The hydraulic 

conductivity dependence explains why the tension difference curve (Figure 23A) changed 

gradually with depth while the moisture difference curve (Figure 23B) decreased abruptly 

between 25 and 40 cm. The plowed layer consisted of loosely packed coarse silt loam 

3 3 

with a bulk density increasing with depth from 1 .15 g/cm to 1 .37 g/cm , while the lay“ 
er between 40cm and 55cm had finer texture (Figure 25) and higher bulk density (1 .47 g/cm 
to 1 .37 g/cm ). These differences, probably due to either management practices or 
local depositional pattern of the floodplain sediments, may have been the cause of 
different conductivity values. Hydraulic conductivity data are not available here; however, 
the differences in hydraulic properties are also reflected in moisture characteristic curves 
which are shown (Figure 25) for depths of 10 cm, 25 cm, 40 cm, and 55 cm. For the 
10 and 25 cm corves, a small change in moisture tension between 100 cm and 400 cm 
results in a targe change of soil moisture content; the same tension change will cause a 
much smaller moisture shift at 40 or 55 cm. The above considerations suggest that while 
tension gradients similar to those in Figure 23Amay develop in more than one soil in a 
given area (due to their partial response to atmospheric evaporative demand), moisture 
contents in these soils can be quite different. The same is true for different depths in 
one profile. For exan^le, tension values on 9/27/73 were identical at 40 cm and 
55 cm (Figure 23A) but the moisture content at 40 cm was 0.047 cm /cm higher than 
that at 55 cm (Figure 23B). Similarly, whereas the moisture tension decreased almost 
uniformly with depth below 55 cm, moisture contents decreased until 100 cm depth and 
then increased. In other words, moisture profiles in the soil may be expected to be 
nonuniform, particularly when the soil heterogeneity existing under natural conditions 
is token into consideration. 
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Rainfall, cm 



Soil Water Tension, cm 


Soil; Eudora Silt Loam 

Depth Particle Size Distribution 1%) 

(cm) Sand Silt Cl 


4 

60.5 

11.1 

4 

59.5 

11.1 

8 

73.4 

23.8 

4 

75.6 

22.1 


Soil Water Content, cmVcm^ 

Figure 25. Moisture characteristic of Eudora silt loam at Four depths. 




A further lllustraMon of fhe differences in moisture dynamics between the upper 
and lower zone is given in Figure 26 which shows the total amount of water in the 
0-32.5 cm layer and 32,5"152.5 cm layer as a function of time . The values were 
computed from measured moisture contents and the constant layer thickness of 15 cm; 
moisture content near the surface was estimoted from the moisture gradient between 10 
and 25 cm. The upper zone was depleted faster during rainless periods (e.g., ^12/73 
to 8 / 30 / 73 ) than the lower zone. However, larger fluctuations in the upper zone 
water storage were caused by rainfall as is evident by comparing the two curves (Figure 
26) for the period between ^31/73 and 9/22/73, During this period, the generally 
low rainfall amounts (Figure 24) did not Infiltrate into the lower zone. The maximum 
range of water storage in the upper zone during the V 12/73 to 9/22^73 period was 
2.085 cm which represented about 25% of the mean storage value; equivalent figures 
for the lower zone ore 0.960 cm and 3%. 

From data presented in Figure 23 and 26 it is evident that soil moisture in the 
upper zone changed more often and reached more extreme values than that in the lower 
zone; in addition, a sharp break in intensity of these changes occurred between 25 and 
40 cm below the surface. Precipitation amounts were identified os one reason. Figure 
27 shows the depth of penetration as a function of the piecipitation amount; the least 
square second-order polynomiol curve indicates that a daily rainfall of 1.27 cm or 
less would not penetrate below 30 cm. During the 106-day long experimental period, 

33 days with precipitation above 0.0254 cm were recorded in Lawrence, Kansas, with 
a total of 48.92 cnv this yields an average daily rainfall of 1 .48 cm, and a correspond- 
ing depth of penetration of 34.8 cm. Consequently, even during this unusually wet 
period the precipitotion amounts were not large enough, on the average, to significantly 
affect moisture contents in the lower zone . 

The second reason behind differences in soil moisture dynamics is evaporative 
loss of water to the atmosphere. The sharp break in the moisture difference curve 
(Figure 23B) indicates that the evaporation loss was small at the depth of 40 cm . A 
further approximation of the lower zone contribution to actual evaporation was obtained 
by the following procedure based on the water balance technique: 

(i) calculation, fa' a rainless period 8/14/73 through 8^30/73, of daily 
drainage out of the soil profile (i.e,, below the lower zone) from 
Equation (6); 

(ii) calculation of total daily loss of water from the lower zone (32.5 cm 
to 152.5 cm) only. 
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Total Soil Water Content, cm 



Depth of Penetration D, 



Results of (i) and (li) are plotted for corresponding days in Figure 28, If the 
evaporation contribution from the lower zone were zero, the points should fall on 
the 1:1 line. However, most of the points are located above the line, suggesting 
that some water from the lower zone moved toward the surface and was lost through 
evaporation; the scatter of points precluded a reliable estimation of the proportion 
of actual evaporation. It should be noted that some points are not shown on Figure 
28 because fbeir coordinates were too large for the scale shown . 

It may be argued that the smalt evaporation contribution from the lower zone 
was caused by the fine~textured layer at depths 40-55 cm which restricted water 
movement. However, Jackson's (1973) data also show that for the first 16 days 
(March, 1971) and 7 days (July, 1970) after irrigation with 10 cm of water, an average 
of 55% of evaporated water originated in the top 9 cm of the Adelanto loam (deep, 
fairly uniform soil). Furthermore, layers with reduced capability to conduct water are 
common in agricultural fields (Saver et al ,, 1972), Trouse and Saver (1965) found 
that almost every agricultural implement created compacted layers under moist condi- 
tions which decrease the rate of water movement through the soil. Consequently, the 
conclusions arrived at for Eudora silt loom appear valid for a wider class of soils. 

2,3,2 Moisture Changes in the Lower Zone 

Since evaporation loss from the lower zone is relatively small, the only 
components which enter the water balance equation (Equation (7)) fw the lower 
zone, in addition to Initial water storage, are precipitation which infiltrates through 
the upper zone, and drainage which escapes below the terminal depth of the soil profile. 
Precipitation is the main source of variations of the moisture stotus in the lower zone; 
consequently, the depth-time changes in soil water content following rainfall were 
of main interest here . 

During a rainless period, moisture losses in the lower zone occur but are smoill; 
these losses are mainly due to drainage below the terminol depth. Figure 29 shows three 
moisture profiles for an almost rainless period between 8/20/73 and S/30/73, Moisture 
decrease at depths below 40 cm was less than 0,01 cmVcm^ for this 1 1-day period, 
while the upper zone was loosing water continuously. Precipitation on 8/31/73 
(3,27 cm) and 9/2/73 (0,78 cm) raised moisture content up to a 4(hcm depth by 


Moisture Loss from 32 . 5 to 132. 5 cm Zone, cm 


▼ 



Figure 28. Lower zone contribution to the actual evoporotion. 


Soil Water Content Change, 

cmVcm^/Period Soil Water Content, cm%m^ 



Figure 29. Soil water content (right-hand side) and soil water content changes (left-hand side) 
before and after rainfalls of 3.27 cm (8/31/73) and 0.78 cm (9/2/73). 


9 / 2/73 bu^ lower depl'hs remained unchanged. Between 9 / 2/73 and 9 / 6 / 73 , moisture 
decreased rapidly above 25 cm, slowly between 25 and 40 cm, and increosed some" 
what at greater depths as a result of water redistribution from the upper layers; this 
increase was largest at 55 cm and smaller below, but it did not exceed 0.010 cmV 
cm . The left part of Figure 29 shows moisture changes at individual depths between 
the dates mentioned. It is apparent that most of the moisture dynomics was confined 
to the upper zone and a few days after the rainfall . 

Whereas the small rainfall reached only to a depth of 40 cm, a large precipi" 
tation on 10/11/73 (8.61 cm) penetrated through and increased moisture at all depths 
(Figure 30). However, more than half of this newly acquired moisture was lost through 
drainage during the next three days (left part of Figure 30). On 10/20/73, lower 
zone moistures were at the levels of } 0 / 6 / 73 , i.e., prior to the large precipitation. 

The rate of moisture decrease further declined as is apparent from curves for the periods 
10/16 to 10/20 and 10/20 to 10/31 , On 10/31/73, the amount of water in the lower 
zone reached a level typlcol for prolonged rainless periods (Figure29). 

These results thus suggest that after the lower zone moisture content was increased 
through precipitation, water was rapidly redistributed to greater depths or lost out of 
the profile; the process was most intense after the addition of water and less so as the 
time progressed, furthermore, the moisture contents tended to assymptotically approach 
a relatively stable vertical distribution in the profile; once this condition v/os tsoched, 
the lower zone moisture changed only slowly. 

Richards et al. (1956) conducted on experiment on Pachappa sandy loam the 
results of which support the above conclusions. They measured moisture content changes 
in the 0 to 50 cm layer for 60 days after an initial irrigation by 45 cm of water. Curves 
in Figure 31 were generated from equations derived by Richards et al . (1956) . The curves 
show that moisture loss was greatest near the surface (due to evaporation) and shortly 
after irrigation but decreased exponentially with time. After 20 days, the loss at 
depths greater than 30 cm was less than 0,0015 cm /cm /day. The exponentiol 
decay type"relotionship for moisture loss due to drainage was also confirmed by 
Wilcox (1959) and Ogata and Richards (1957); evaporation was excluded in these 
studies. 

The conclusions regarding rapid redistribution rate and establishment of an 
assymptotical vertical moisture distribution are valid only if sufficient amount of 
water is available initially and provided that water movement in the soil is not 
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Depth, cm 


Soil Water Content Change, 

cm%m^/Period Soil Water Content, cm%m^ 

„-0.050 0.000 0.050^ 0.200 a 300 



o 10/6-10/11 
10/11 - 10/14 
•10/14-10/16 
• 10/16 - 10/20 
■ 10/20 -10/3L 




Lcatr* 


rear* 


Eudora Silt Loam 




Date-. 

o 10/6/73 
10/11/73 
• 10/14/73 
▼ 10/16/73 
■ 10/20/73 
□ 10/31/73 


Figure 30. Soil water content (right~hand side) and soil water c<»itent changes (left~hand side) 
before and after a rainfall of 8.61 cm (lQ/11/73). 


Pachappa Sandy Loam 

Depth: 


• 

0-10 cm 

▼ 

10-20 cm 

o 

20 -30 cm 



30-40 cm 


40-50 cm 


Days Since irrigation 


Figure 31 . Soil water content change as a function of depth ond time 
Based on the data of Richards et ol. (1956)* 











appreciably inhibited. In that cose, soil releases all water which it cannot retain. 
However, in many situations, soil is capable of retaining more water than is available. 
An example is water depletion by plants. Figure 32 presents average monthly soil 
moisture storage for four geographic locations (Mather, 1964) computed according to 
the climatic water balance technique of Thornthwaite and Mather (1957); the 
computations assumed soil moisture storage of 30 cm. At three of the four locations, 
the soil water is not replenished during the entire year. If water were supplied here, 
it would not be completely released as drainage until the soil storage capacity is 
replenished, in one cose (Topeka), the rnoistu-e storage is full for four months and In 
this period the conclusions reached previously would apply. Secondly, the data 
referred to in this section were measured in soil profiles with apparently high hydraulic 
conductivity which facilitated rapid movement of water. The rate of response to preci" 
pitation would be slower in the case of soils with low conductivities,e.g. clays. These 
considerations are important because they place restrictions on the algorithms that 
might be designed to estimate the actual moisture storage of the lower zone. 
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Soil Profile Water Storage, cm 


0 Topeka 

39® 04* N 

95° 38' W 

267 

1905 -t/8 

• Austin 

30® 18' N 

97° 42' W 

187 

1905-69 

▼ Bismarck 

46®46'N 

100“ 45' W 

503 

1905 -60 

■ Phoenix 

33® 26' N 

112° or W 

•338 

1905 - 60 



Month 


Figure 32. Mean monthly soil profile water storage at four locations. From 
Thomtwaite and Mather (1964). 




CHAPTER 3. 


SOIL TEMPERATURE REGIME 

Temperalure of a bare soil is a resulf of ini'eracfions of numerous variables, 
both exiernal and int'ernal to the soil profile. Because of Its importance to mass and 
energy ■exchange near the earth's surface, soil temperature and its changes with these 
variables Save been studied extensively by soil physicists, meteorologists, climatologists, 
civil engineers, and others. For example, Willis (1964) listed 1,152 references dealing 
with soil temperature. More recently, attempts have been made to combine the various 
kinds of information into simulation models. If the established concepts are correct, 
computed temperatures should agree with the measured ones. Once the validity of a 
model is established, the model may be used to piedict the effect of changing various 
parameters, environmental conditions, etc., whereua it would be very difficult to 
collect equivalent data experimentally. 

Soil temperature changes with time are described by the diffusion equation, and 
the analytical opproach to soil temperature modeling has been based on simplified so- 
lutions of this equation. Several examples of this type were given by van Wijk (1963) ,vanWi]k 
and de Vries (1963), and van Wijk and Derksen (1963). These solutions hold for cases 
in which soil is homogeneous and boundary conditions can be described analytically. 

Under natural conditions, however, soil thermal properties vary with depth and the 
upper boundary condition (surface temperature) changes may deviate from ancnalytical 
function. For these reasons, numerical modeling has been attempted as computers became 
available; its advantage is that the two above restrictions are no longer necessary. 

Numerical models involving soil temperature have been employed by soil 
physicists (Wierenga et al., 1969; Wierenga and de Wit, 197Q; Hanks et al., 1971; 

Westcot and Wierenga, 1974), and meteorologists or climatologists (Effimova and 
Tsitsenko, 1963; Myrup, 1969; Outcalt, 1972; Foster, 1972; Foster and Fye, 1973; 

Dyer, 1974). The former studies are of limited interest here because they use a 
measured upfjer boundary condition, usually 1 cm or so oelow the surface, in contrast, 
the latter models compute surface soil temperature and then use it as a boundary condi- 
tion for the heat wave propagation below the surface. Most of the above meteorological 
models compute surface temperature from the energy balance equation; differences exist 
in the assumptions made and input parameters . 
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To determine fhe magnitudes of bare soil temperature fluctuations under variois 
conditions, a numerical model of soil temperature was prepared (Appendix A); it 
followed, in its general features, the model developed by W, D. Sellers and described 
by Foster {1972}. 


3.1 SOIL TEMPERATURE MODEL AND ITS PERFORAAANCE 


The energy balance equation for the soil surface can be written as 


= H + LE + G. 


06) 


is net radiation and consists of absorbed solar (shortwave) radiation and absorbed long- 
wave atmospheric emission minus thermal infrared radiation emitted by the soil surface. 

H and LE are sensible and latent heat terms, and G is soil heat flux. All tenns in 
Equation (16) can be positive or negative. If the aerodynamic concept is used to 
compute H and LE, an analytical solution of the diffusion equation to compute G, end 
semiempirical equations to calculate the components of Rj^, soil surface temperature Tg 
can be found using Equation (17); 


T = 


where 


T + 
a 


S -Co-C.-C, 
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= absorbed solar radiaMon, in !y/ min; 

T = air t■empera^ure, in ®C; 

Q * \ 1 A 

O' = Stefan“Bolf'zman constanf, a = 8,17 * 10 ly/min 
^ = infrared emissivily of fhe soil, dimensionless; 
n = cloud cover in fenfhs, dimensionless; 

k = paramefer varying wiih cloud type, given by Sellers (1965, p. 58); 
p = air pressure, in mb; 

Cp= specific heat of air, in cal/g/^K; 

D = transfer coefficient, in err/ min; 

R = gas constant, R = 2870.4 mbcm / °K/g; 

L = latent heat of vaporization, in cal/g; 

a = AE/PE ratio, dimensionless; 

e^^ = saturation vapor pressure at T^, in mb; 

= actual vapor pressure at the height of wind speed measurements; in mb; 
to = angular frequency, in radians per time increment of soil temperature 
calculation; 

C = soil heat capacity, in cal/cm / C; 

^ = soil thermal conductivity, in cal/cnr/sec/°C; 

T = average doy soil surface temperature, in °C; 

^ o 

= soil surface temperature for the previous time increment, in C. 

The derivation of Equation (17) is given in Appendix A. 

Assuming horizontal uniformity and heat transfer by conduction, the temperature 

change with time t at depth z, T(z,t), is described by 

( D*(z) fe-fil ) (18) 

6 t Irz ^ ^ • 6z ' • 

2 

where D'(z) = apparent thermal diffusivity in cm /min. The computer program incorporaf 
ing Equation (17) and the finite difference form of Equation (18) is also described in 
Appendix A . 

Before the mode! can be used to simulate soil temperature variations under 
various conditions, its performance must be tested to ascertain its correspondence to 
reality. Two sets of profile soil temperatures measured over a 24 hour period* were used 

*Courtesy of Dr. S. B. Idso, U, S, Water Conservation Laboratory, Phoenix, Arizona. 



for rhis purpose; fhe experimental procedure used for rheir acquisition Is described 
In Appendix B. Table 10 gives volues of bulk density, field capacity, initial 
moisture content and soil temperature for both sets; it should be noted that bulk 
densities below 10 cm and soil moisture below 30 cm were assumed as no direct meo“ 
surements were available. Field capacity was taken as moisture content at 110 cm 
tension (Jackson, 1973), and the value for 0-2 cm was reduced by 0,9 (see Section 
2.2.1). Additional variables for both test runs ore listed in Table 11. 

The correspondence between computed and measured soil temperatures for the 
9/2Q/7Z set (wet soil) is shown in Figure 33 and 34. Calculated surface temperatures 
followed the measured values fairly closely, although a systematic shift occurred: the 
computed surface temperatures increased and decreased more slowly than the measured 
ones (Figure 33). The difference was greatest at 1600 hours (3.5®C) but decreased to 
small values during night hours. The diurnal variation of computed and measured 
temperatures at 5 cm, 10 cm, and 20 cm (Figure 34) also indicates reasonable 
predictive occuracy of the model. At all depths, however, the predicted diurnal 
temperature amplitude was greater than the measured one. 

Comparison of computed and measured soil temperatures for dry soil (](y 02/73) 
revealed similar results (Figures 33, 35). Surface temperature was predicted accurately 
until the maximum value at 1300 hours, but It decreased more slowly than measured 
temperatures between 1300 and 2000. Subsurface predicted temperatures also followed 
the meosured values until approximately 1300, after which the difference between the 
two increased. 

Results shown in Figures 33 through 35 suggest that the accuracy of soil temperature 
prediction differed somewhat for surfoce and subsurface temperatures. The lag in surface 
temperatures of the wet soil (Figure 33) during the morning hours could be due to the 
latent heat; this possibility is enhanced by the differences being negligible for the 
dry soil case (Figure 33). Similarly, an error in predicting the latent heat term could 
be responsible for the lag of computed surface temperatures In the afternoon; note that 
on 10/02/73 for example, the moisture loss proceeded faster than the model predicted, 
thereby leaving less heat Fot warming the soil (depth 2 cm, time 1300 hours to 1800 hours. 
Figure 15), In general, however, surface soil temperatures were predicted with good 
accuracy as the obove mentioned differences were not large , In the case of subsurface 
temperatures, diumol temperature amplitude was consistently larger. Indicating that 
the computed heat conduction into the soil was too high. The apparent thermal 
diffusivity employed in computing temperature changes from Equation (18) was 
determined using thermal conductivity calculated from de Vries' (1963) model, and a 
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Table 10. Bulk densify, field capdclfy, initial soil moisture, and initial soil 
temperature profiles for two soil temperature data sets used for a 
soil temperature simulation model testing. 


Depth 

(cm) 


10/02/73 


9/20/73 & 10/02/73 


Bulk 

density 


(g/cm^) (cmVcm^) (cm^cm^ (°C) j(cmVcm^j (°C) 


Moisture 

Tempera- 


ture 

(cmVcm^ 

(°C) 









Table 10. Bulk density, field capacity, initial soil moisture, and initial soil 
( Continued) profiles for two soil temperature data sets used for a 

soil temperature simulation model testing . 


Depth 

(cm) 


9/20/73 & 10/02/73 


Bulk 

density 


(g/cm^) 



0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 

0.305 


9/20/73 

IQ/02/73 

Moisture 

Tempera- 

ture* 

Tempera- 

ture** 

Moisture 

Tempera- 

ture 


(°C) 

("*C) 

(cmVcm^ 

(°C) 

0.223 

29.7 

29.5 

0.240 

27.9 

0.220 

29.9 

29.6 

0.240 

27.9 

0.220 

30.0 

29.8 

0.240 

27.9 

0.220 

30.1 

29.9 ■ 

0.240 

27.9 

0,220 

30.2 

30.1 

■ 0.240 

27.9 

0.220 

30.3 

30.2 

0.240 

28.0 

0.220 

30.5 

30.4 

0.240 

28.0 

0.220 

30.6 

30.6 

0.240 

26.1 

0.220 

30.7 

30.7 

0.240 

28.1 

0.220 

30.8 

30.9 

0.240 

28.1 

0.220 

30.9 

31.1 

0.240 

28.1 

0.220 ' 

31.0 

31.3 

0.240 ' 

28.1 

0.220 

31.1 

31.5 

. 0.240 

28.1 


Note: 


Measured vales. 


it 

Values computed by starting with a constant soil temperature with depth, and 
then computing successive iterations until the temperatures at depth 0 cm and 
70 cm at 2400 hrs, differed by less than O.rC for two adjacent diurnal runs. 
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Table 1 1 . OperaMonal parameters for the soil temperature simulation model . 


Data Set 


Hypothetical 


Parameter 



‘VI — 


9/20/73 

10/V73 

HMP* 

LMP* 

PERL 

10 min. 

iO min. 

10 min. 

10 min. 

DT 

2.5 min. 

2.5 min. 

2.5 min. 

2.5 min 

Number o. layers 

35 

35 

35 

35 

Anemometer height 

60 cm 

60 cm 

240 cm 

240 cm 

Roughness length 

0.003 cm 

0.003 cm 

0.5 cm 

0.5 cm 

Previous day average 

27®C 

25°C 

26,2°C 

28.6°C 

surface temperature 





Previous hour surface 

21®C 

18°C 



temperature 





Emissivity of the soil 

0.90 

0.90 

! 

0.90 

0.90 




HMP (LMP) = high (low) moisture content profile 





e 
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Soil Temperature. ®C 


Adeianto Loam 
9 / 20/73 



Figure 34. Comparison of computed and measured soil temperatures at three depths/ 
wet soil. 



Time of Day, Hours , ] 



Figure 35. Comparison of computed and measured soil temperatures at three depths, 
dry soil. 





simple formula was used for esfimaring heot- capacity (see Appendix A). Both of 
these methods have previously been found to yield satisfactory results (Hanks et al 
1971; Wierenga et al,, 1969; de Vries, 1963). The thermal diffusivlties generated in 
this manner for the dry soil case ranged from 0.353 cmVmln. at 2 cm depth to 0.410 
cm /min at 10 cm. Considering the relatively high bulk densities of Adplanto loam 
(Table 10), these values are consistent with those determined by Wierenga et al. 
(1969) as well as those used by Hanks et al. (1971) for similar conditions. To estimate 
what the diffusivity should be in order for the computed temperatures to follow the 
measured ones, an "effective” thermal diffusivity tx was computed from the dry soil 
data using an amplitude equation (Wierenga et al., 1969); 

“ = < MA/ ?7r) ' 

where A^ (A 2 ) is the amplitude of the diurnal temperature wave at depth z^(z 2 ). The 
following results were obtained; 


Z] (cm) Q 

1 

2 

3 

4 5 10 

0 

0 

^2 (c»«) 1 

2 

3 

4 

5 10 20 

5 

10 

A^(®C) 39.5 

29.1 

23.2 

19.3 

17.1 15.2 10.5 

39.5 

39,5 

£H(cmVmin) 0.0234 

0.0425 0.0644 

0.1490 

0.1573 0.3986 0.2738 

0.0598 

0.1243 


The computed diffusivlties in the 0*5 cm layer are unusually low. For example, assuming 
heat capacity of 0,35 cal/cmV^C, a diffusivity of 0.0598 cm^min (depth 0 to 5 cm) 
means that the apparent thermal conductivity would have to be 3,49 x 10 ^ cal/cn/sec/ 
®C; this value is 60 times smaller than the thermal conductivity of quartz at 10°C given 
by de Vries (1963), Some of the discrepancy is probably due to the fact that Equation 
(19) assumes constant diffusivity with depth. Considering the dry soil vertical moisture 
distribution (from approximately 0.02 cm /cm at 0 cm to 0.17 cm /cm at 5 cm), 
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this assumption was clearly violated. Nevertheless, the rapid decrease in the 
temperature wave amplitudes with depth indicates that heat penetration into deeper 
layers was slrongty attenuated. If heat transfer is assumed to have occurred by 
conduction only (as described by Equation (18)), theii the above observations imply 
that the computed thermal diffusivities must be reduced in order to achieve satisfactory 
correspondence between computed and measured temperature values. Consequently, 
the thermal conductivities computed from de Vries* (1963) model were multiplied by 0.3 
in the dry soil case, while no correction was applied to wet soil computations. Figure 
34 shows that the reduction was adequate for the period of low heat load only. This 
suggests that other processes may not have been sufficiently accounted for, e.g. latent 
heat loss. 

In spite of these discrepancies, the prediction of the soil temperatures was 
reasonably satisfactory, the computed values for most time/depth/moisture combinations 
being within 1 .5°C of the measured values. In the following sections, these discrepancies 
will be of minor importance because relative changes in soil temperature will be of 
primary interest. 


3.2 DIURNAL SOIL TEMPERATURE CHANGES 

The soil temperature simulation model described and tested in the previous section 
can now be used to estimate the effect of importont environmental variables on soil 
temperature. Since changes in soil temperatures ranging from very small to very large 
con be caused by numerous parameters, some restrictions must be imposed to make the 
analysis manageable. Consequently, the problem here is limited to estimating pr<^able 
"average maximum” ''oriations of bare soil temperatures for mid-latitudes. 


From t'he many parameters related to diurnal soil temperature variations, the 
following are of major importance: moisture content, time of day, zenith ongle at 
solar noon, and depth. Clear sky is assumed because then extreme temperatures can 
be reached. To render results of this analysis more representative, hypothetical 
experimental conditions were defined as follows: 

(i) The hypothetical bare field was located near Topeka, Kansas. 

The site had three surface configurations: 

slope 0°; 

slope 30°; aspect +90°, i ,e . facing West; 
slope 30°; aspect ”90°, i.e. facing East. 

Roughness length was taken as 0,5 cm, a value intermediate between 
0.03 cm given by Pasquill (1962) for level desert and 2.1 cm for fallow 
field (van Wijk and Borghorst, 1963). As will become apparent later, 
choice of another roughness length would result in somewhat different 
temperatures . 

(ii) The experiment took piece around July 20 when solar declination was 
20.85°N. 

(iii) Hourly values of screen height air temperature, wind speed, and 
relative humidity were obtained by averaging observations taken at 
3“hour intervals during four clear sky days in July, 1973, at the 
Topeka Weather Station (latitude 39°04’N, longitude 95°38'W, 
altitude 267 m). Solar radiation was calculated (Appendix A) assuming 
the coefficient of atmospheric transparency in Kostrow’s equation 

TC = 0.314 for the optical mass OM = 1 .0. Under these conditions, 
total solar radiation was equal to 1 .414 cal/cm /min for a horizontal 
surface at the solar noon. The same meteorological conditions were 
assumed for both wet and dry soil . 

(iv) Bulk density was assumed to be 1.25 g/cm throughout the soil profile. 

The soil mineral fraction was assumed to consist of 30% quartz and 
70% other minerals. Thermal conductivity computed by the method 

of de Vries (1963) was not reduced for the dry soil case, primarily because 

the computations yielded thermal diffusivities comparable to those of Hanks 

et al . (1971) andWierenga et al.(1969) for equivalent conditions. Field 

3 3 

capacity was assumed equal to 0.300 cm /cm except for the top 2 cm where 
it was reduced by a factor of 0.9 (see section 2.2.1). For the high moisture 


84 


con^en^ profile (HMP), soil was considered fo inifiolly confoin 
0.300 cm^crn^ i'hroughout the profile. The dry profile (IMP) 
was derived by applying the layered water balance model to the 
soil initially at field capacity; the LMP was taken as the moisture 
profile afte«‘ 20 days of drying without precipitation during which 
potential evaporation was a constant of 0,575 cir/day. 

(v) Initial soil profile temperatures were derived by means of the soil 
temperature simulation model, separately for HMP and LMP. Using 
the above conditions and an initially constant soil temperature, 
iterrations were computed until temperatures at 2400 hours changed 
less than 0,1°C for depths 0 cm and 70 cm between successive 
iterrations. The adequacy of this procedure was verified on the 
9/2Q/73 data set (see Table 10). 

Diurnal soil temperatures discussed below were computed during six different 
runs. The results are identified by three parameters, namely moisture condition (H = 
high, L « low), slope (in degrees), and aspect (in degrees); for example, H/30/“90 
refers to temperatures of soil with high moisture content, surface slope 30°, and surface 
facing East, 

Figure 36 shows diurnal changes in surface soil temperature for the two moisture 
contents (H/O/Oand L/0/0) and roughness length 0.5 cm. Before dawn, surface 
temperature of the wet soil was slightly higher, apparently due to higher heat 
conduction from deeper layers (Figure 39, 40). Following sunrise, dry soil surface 
temperature increased much faster than thot of the wet soil as a result of lower latent 
heat expenditure. Maximum surface temperatures were reached approximately 80 
minutes ofter solar noon in both coses which agrees with the observations of Idso et al . 
(1975). The surface temperature wave amplitude was doubled for dry soil compared to 
wet soil . Soil moisture thus exerted a major effect on the daytime surfoce temperature. 
However, Figure 36 suggests that this amplitude is also a function of surface roughness. 
The pairs of curves were computed for identical conditions except for the roughness 
length which was 0.50 cm and 0.01 cm, respectively. As a consequence of smaller 
roughness, the surface temperature amplitudes increased by a factor of 1 .5 for both 
wet and dry soil. This difference was caused by lower turbulent transfer coefficient 
D in the smooth case which decreased latent heat transfer. For example, for the wet 
soil, latent heat decreased from 516,6 ly/day (rough) to 376,0 ly/day (smooth). 

Surface temperature thus increased, in spite of an increase in sensible heat exchange 
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(from "48.0 \y/day for rough surface to 26,9 ly/day for smooth surface). It should be 
noted that the surface temperature computed by the mode! represents an "averoge” value 
over the surface soil elements. Soil moisture content also affected soil temperatures at 
greater depths; this is evident frorri Figure 3!^ and 38b which show temper ati res at 10 cm. 

With Increasir^ surface temperature amplitude, the amplitudes ot greater depths also 
increased. This occurred because the thermal diffusivity was not appreciably reduced | 

for the IMP (both thermal conductivity and thermal capacity decreased, resulting in a i 

small diffusivity change). 

The effect of slope and aspect of soil surface is shown in Figure 37 for the dry 
soil. As expected, surface temperature was highest for l/30/“90during morning hours 
and for L/30/90 in the afternoon; this was determined by the actual solar zenith angles 
for individual surfaces. The daily maximum surface temperature increased somewhat os 
the aspect changed from east to west, but the increase was small for LMP (Figure 37) j 

and negligible for HMP (Figure 38). These data suggest that the surface temperature 
may vary considerably example, surface tempetahjres L/30/90 and L/30/“90 | 

differed by 13.3 °C at 0900 hours. The differences were somewhat lower in the after- | 

noon. The slope/aspect effect was smaller for the HMP (Figure 38) because of the lower 
surface temperature amplitude but the trends were similar. In both cases, the differences 
in surface tempetatures were propagated to greater depths (Figure 37,38). 

The changing depth/time pattern is illustrated in Figure 39 (L/O/0) and Figure 
40 (H/C/O). Between midnight and dawn, soil temperature increared almost linearly to 
20-30 cm and remained approximately constant at greater depths. Following sunrise, 
surface temperature increased first and the heat wave began to penetrate to deeper layers. 

After surface temperature started to decrease, the temperature profile had both concave 
and convex portions, causing heat transfer to proceed in two directions. This condition 
was well developed by 1800 hours, and from then on soil temperature decreased. The 
time/depth patterns for dry (Figure 39) and wet (Figure 40) soil were similar, the main 
difference being temperature amplitudes at different depths. The amplitudes shown 
in Figure 41 permit two observations: (i) on the diurnal basis, soil tempetature changed 
relatively little at depths below 30 cm; and (ii) the amplitudes were doubled in the cose 
of dry soil . 

Data for the above soil temperature simulation were chosen so as to represent 
"average maximum" changes in bare soil temperature to be expected at mid-latitudes. 

The good performance of the model in matching measured data suggests that the simulated 
results are representative for the chosen environment. Considering all factors that affect 
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Figure 38. Temperature of o wet soil at the surface (A) and at 10 cm depth (B) as affected 
by the surface configuration. 
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Figure 39, Temperature of o dry soil as a function of depth and time of day 
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surface temperature and magnitudes of the changes, it is apparent that it would be 
quite difficult and impractical to predict bare soil temperature for small parts of an 
extensive area. Such a prediction would be further complicated by variable cloud 
cover. However, the accuracy to which soil temperature should be krown for micro- 
wave remote sensing of soil moisture is not obvious at this point; therefore, this 
question will be addressed again ofter the soil temperature date are analyzed in terms of 
the microwave signal In Chapter 4. 


CHAPTER 4. 







SOIL MOISTURE AND TEMPERATURE REGIMES AS A FACTOR 
IN MICROWAVE REMOTE SENSING OF SOIL WATER CONTENT 

In ^he previous secfions, the attention was focused on soil teir *ature and 
moisture regimes per se, without studying their implications regarding microwave 
remote sensing of soil moisture. Mathematical models were employed to illustrate 
major features of these regimes as well as to generate hypothetical sets of data. 

These models allowed examination of some features of the moisture and temperature 
regimes under relatively extreme conditions. In addition, actual meosuremenH 
collected during various field moisture experiments were analyzed to examine 
particular aspects of the moisture regime. Results arrived at in the previous chapters 
con now be evaluated from the viewpoint of the microwave method of remote sensing 
of soil moisture. A microwave sensor designed for remote sensing of soil moisture 
would operate in active (radar) or passive (radiometer) mode and at specified frequency, 
incidence angle, and polarization. To render the analysis manageable, prim jry 
attention is given here to three frequencies (1 .4 GHz, 4.0 GHz, 10,0 GHz; 1 GHz = 
10^ cycles/second), normal incidence (incidence angle 0®),and a passive system. The 
frequencies were chcsen for two reasons. First, dielect, ic constant measurements 
were available at 1,3 GHz, 4,0 GHz, and 10.0 GHz (Qhlar and Ulaby, 1974). 
Secondly, the 1 .4 to 10,0 GHz ronge encompasses most of the recent experimental 
studies (Ulaby, 1974; Ulaby et al,, 1974; Batiivaio and Cihlor, 1975; and Newton 
et ol., 1974) and thus a microwave sensor for soil moisture sensing would likely 
cover one or more of the intermediate frequencies. The discussion is limited to 
norma) incidence in order to simplify the analysis of polarization and angular effects. 

As shown in Chapter I, the ptissive (radiometric) signal is directly influenced by both 
temperature and moisture of the soil; therefore, the passive mode was chosen here to 
illustrate the combined effects of the previously discussed soil moisture and te* iperature 
regimes. In addition, a limited comparison is made between power reflection 
coefficient (active case) and emissivity (passive cose) values computed independently 
for identicol soil moisture profiles, ond the effect of soil temperatures on soil 
dielectric constant (and therefore directly on power reflection coefficient and 
emissivity) is illustrated for relatively extreme conditions. More extensive analysis 
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of the active case was limited by financial constraints as the power reflection 
coefficient computations were rather costly. It should be emphasized that 
although the above choice of mode or system parameters was made with the 
objective of accurately portraying the effects of soil moisture and temperature 
regimes on the microwave signal, he choice was arbitrary and not intended to 
imply optimum sensor design criteria. 

The second problem addressed in this chapter was the possibilit)' of estimat~ 
ing moisture content within the soil profile and over large geographical areas. The 
algorithm proposed in section 4.2 is based on the assumption that a microwave 
remote sensor is available which yields near^surfbce soil moisture information on a 
periodical basis. Again, the optimum sensor characteristics have not been dealt 
with here, as an entirely different approach to the latter problem would be required. 


4.1 EFFECT OF SOIL MOISTURE AND TEMPERATURE REGIMES ON THE 
BRIGHTNESS TEMPERATURE 


According to Equation (2), the radiometric temperature T^p is a function of 
soil emissivity e at a given frequency, physical soil temperature T^, and of the 
radiation T emitted by the upper hemisphere; T consists of contributions by the 
gaioxy (Tg^j^ ) and of the atmosphere (T^j^jj), respectively. Since the microwave 
energy measured by the radiometer is emitted at various depths in the soil prior to 
emerging above the surface, and (Equation (2)) may be considered "effective" 
quantities describing the total soil effect on the microwave radiation. If T and e 
are constant with depth, then they will be equal to these effective quantities and the 
brightness temperature T^^ can be readily determined. However, both moisture 
(and therefore emissivity) and terr^erature vary with depth. TTierefore, in order to 
determine the effective T^ and e values in this case, some kind of weighting as a 
function of depth must be employed. 

Burke and Paris (i974) described a model based on radiation transfer theory 
which predicts brightness temperature over a bare soil. This model assumes that: 

(i) soil consists of a number of layers with smooth parollel interfaces 
between adjacent layers; 

(ii) radiation is inedierent; 
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(iii) there is no attenuation or emission between the surface and 
the sensor; 

(iv) the sky brightness is isotropic; and 

(v) dielectric properties are constant across any given layer 
of the soil . 

th 

If the soil consists of N layers (N layer has an infinite thickness), then brightness 
temperature immediately above the surface measured at a given frequency, polarization q 
and look angle 6 is equal to 


N 

^bq ~ ^pi (e)Az.))( 1 + R {0)exp(-Y.(e)Az.)) 

i==i 


I I 

• 1 1 ( ^ “ R : (e ))exp (- 2 Y , (e)Az , ) , 

j=l £,=2 ^ 


( 20 ) 


Y5(e) 


2 0 } a .(e) 
21 ^ ^ 


C 


a , 

21 



f 


B . = 

zi 



sin^( 0 )) r 
1 + 



where 

Az. 

R .(0) 
qj 

c 

0) 

k* 

k" 


= thickness of the i*^^ layer, in cm; 

= power reflection coefficient for a given polarization 
q, layer j , and look angle 6, dimensionless; 

= velocity of light, in crr/sec; 

= angular frequency, Hz; 

= reol part of the relative dielectric constant, dimensionless; 

= Imaginary part of the relative dielectric constant, dimensionless. 
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A computer program written to compute EquoHon (20) was provided for the present 
study.* The program was modified to fit data generated in Chapter 3 and expanded 
so that the following parameters could be calculated for each moisture/temperature 
profile combination; 



( 0 ) 




( 21 ) 


A 6 




(e) 


( 22 ) 


e 

q 


( 0 ) = 




( 0 ) , 


(23) 


m 


eb 



■"i «» 

Tb,(0> 


(24) 


N 



"■i ^ ^ qj 

( 0 ) ' 


( 25 ) 


*Courtesy of Dr. W. J . Burke, National Research Council, Johnson Space 
Center, Houston, Texas* 
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where 




= brighfness femperal'ure for pola^1za^io^ q immedial'ely above 




qj 


'■q 

tn 


eb 


ea 


m . 
J 


the surface, in K; 


^bqj 


?th 


= brighiness ieniperal'ure contribufion of the j ” layer for 

polarization q when it emeiges above the soil surface; 
til 

= j layer contribution to the total emissivity dimensionless; 
= total {i.e. "effective") emissivity of the soil, dimensionless; 

= sum of moisture contents of individual layers weighted by the 

contribution of individual layers to the total brightness 

.33 

temperature, in cm /cm ; 


m = sum of moisture contents of individual layers weighted by the 


3 3 

contribution of each layer to the total emissivity, in cm /cm ; 

th 3 3 

= volumetric water content of the j layer, in cm /cm . 


Reasons for defining two different measures of equivalent moisture will become obvious 
later. For the present, note that in the case of constant temperature and moisture 
profiles with depth, m^j^ = ”'ea* 

Equations (20) through (25) were applied to soil moisture and temperature profiles 
computed in Chapter 3 for the following coihbinations: 
time; 0000, 0600, 1200, 1800, 2400 hours; 
moisture; high (H), low (L) water content; 
surface configuration (slope/aspect): 0/0, 30/90, 30/“90; 
frequency: 1 .4 GHz, 4.0 GHz, 10.0 GHz. 


Throughout section 4,1, it should be kept in mind that the moisture/temperature 
profile pairs were generated assuming a clear~sky day, around July 20, and for 
mid~latitudes . The above combinations therefore illustrate the diurnal (0000 to 
2400 hrs) changes in moisture and temperature to be expected under these condi** 
tions for initially moist and initially dry soils, respectively. Tables 13 and 14 
contain moisture and temperature profiles for H/O/0 and l/O/O as the representative 
combinations. Dielectric constont values k* and k" were computed from polynomials 
(Table 12) based on data for loam provided by Cihlar and Ulaby (1974). The 1 .3 GHz 
data given by Cihlar and Ulaby (1974) were considered representative of 1 .4 GHz. 
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Table 12. Polynomial coefficienfs for computing the real and imaginary parts 
of the relative dielectric constant of a moist soil .* 


Frequency 

(GHz) 


Coefficient 


Red (k') 3.611374 

Imaginary (k") -0.010487 

Real 3.18690 

Imaginary 0.112123 

Real 2.38487 

Imaginary | 0.23435 


7.45875 

7.11036 

-1.03203 

6.34419 

16.3802 

2.9045 


137.3191 

-12.07477 

163.3028 

18.2682 

83.2034 

63,4922 


me general equotion is 


k' or k" = Oq + a^m + 


3 3 

where m = volumetric soil water content in cm /cm . 









TemperaJore dependence of k* and k" was nof faken info account since most of the 
experimental k' and k" data were obtained at 20°C. Values of Az. = 2.0 cm were 
used for 1.4 GHz and 4,0 GHz while computations at 10.0 GHz were made with 
Az. = 1,0 cm to achieve better depth resolution* Although calculations were made 
for look angles between 0° and 90®, only 0® results are discussed because here the 
effect of soil temperature and moisture on brightness temperature is most explicit. 

At this angle, brightness temperatures for horizontal °tid vertical 
polarizations are equal and will be subsequently referred to as T^(0). Furthermore, 

€ (0) = e (0) = € , and Ae .(0) = Ae. in the following discussion. 

q qj j 

Figure 42 shows diurnal Tj^(O) changes for H/O/O and l/Q/0 at three 
frequencies. The soil initially at field capacity (Figure 42A) exhibited rapid 
rise in the brightness temperature throughout the daytime but changes during the 
night were small . Since soil temperatures at 0000 and 2400 hrs. were very similar at 
all depths (Table 14), the rise in T^j(O) must be attributed to moisture loss due to 
evaporation. The loss was smoll during nighttime (Table 13) which exploins the 
low brightness temperature change. This statement is further supported by the 1/0/0 
data (Figure 42B) where the evaporation was lower and therefore Tj^(0) values at 0000 
and 2400 were much closer. The brightness temperature dependence on moisture 
content implies that between 1800 and 2400, T|^(0) should decrease more than Figure 42A 
indicates. This is because the layered evaporation model used for computing soil 
evaporation losses (Section 3.1) did not account for near-surface recharge due to 
capillary action during late afternoon and evening hours. According to data in Figure 
15, near-surface moisture content at 2400 equalled approximately that at 1400. 
Consequently, if the recharge were accounted for, maximum T^(0) values for wet soil 
(Figure 42A) would be somewhat lower and T|^(0) would probably decrease to about 
215®Kat2400. This Is also supported by Figure 42B which shows that when moisture 
content changed little during the 24 hr. period, the difference between brightness tem- 
peratures.was also small. The three surface configurations (C/0, 30/90, 30/-90) affected 
brightness temperatures mainly during daytime . The maximum differences occurred at 
1200 and were equal to 10.4°C (1 .4 GHzl, 11.1 ®C (4.0 GHz) and 9.3®C (10.0 GHz), 
respectively. It should be noted that the maximum diurnal range of brightness temperatures 
(ij.g,, 38. 9^ at 4,0 GHz, Figure 42A) was mainly due to high evaporation loss which in 
turn was affected by the roughness length used during temperature regime simulotions 
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Figure 42. Diurnal changes of a bare soil brightness temperature for a wet (A) 
and a dry (B) soil at three frequencies and three surface 
configurations. 
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Toble 13. Moisture confenli for fhe hypofheMcal profiles of various 
depfhs and fimes of day. 


MOISTURE PROFILE 


Depth 

in 

cm 

H/0/0 

n 

-- • 

L/O/O 

3 3 

Moisture Content at Time (hours), in cm /cm 

00 

06 

12 

18 

24 

1 

00 

06 

12 

18 

24 

1 

0.300 

0.296 

0.249 

0.188 

0.181 


0.042 

0.042 

0.040 

0.035 

0.035 

3 

0.300 

0.298 

0.267 

0.231 

0.226 


0.107 

0,107 

0.103 

0.095 

0.094 

5 

0.300 

0.298 

0.280 

0.257 

0.254 


0.156 

0.156 

0.152 

0.147 

0.146 

7 

0.300 

0.299 

0.284 

0.264 

0.261 


0.172 

0.172 

0.169 

0.164 

0.163 

9 

0.300 

0.299 

0.286 

0.268 

0.265 


0.181 

0.181 

0.179 

0.174 

0.173 

11 

0.300 

0.299 

0.287 

0.270 

0.268 


0.188 

0.188 

0.186 

0.181 

0.181 

13 

0.300 

0.299 

0.289 

0.272 

0.270 


0.191 

0.191 

0.189 

0.185 

0.184 

15 

0.300 

0.299 

0.290 

0.275 

0.273 


0.195 

0.195 

0.193 

0.189 

0.188 

17 

0.300 

0.299 

0.291 

0.277 

0.275 


0.203 

0.202 

0.200 

0.196 

0.195 

19 

0.300 

0.299 

0.292 

0.280 

0.278 


0.211 

0.210 

0.208 

0.204 

0.203 

21 

0.300 

0.299 

0.293 

0.283 

0.281 


0.215 

0.215 

0.213 

0.209 

0.209 

23 

0.300 

0.300 

1 

0.295 

0.286 

0.284 


0.222 

0.221 

0.219 

0.216 

0.215 

25 

0.300 


0.296 

0.289 

0.288 


0.230 

0.230 

0.228 

0.224 

0.224 

27 

0.300 


0.297 

0.292 

0.291 


0.240 

0.241 

0.239 

0.235 

0.235 

29 

0.300 

0.300 

0.298 

0.295 

0.295 


0.257 

0.257 

0.255 

0.252 

0.252 

31+ 

0.300 

0.300 

0.299 

0.299 

0.298 


0.278 

0.277 

0.277 

0.276 

0.275 
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(0.5 cm, section 3,2). Whether or not such a T|^(0) range would occur under field 
conditions would to a considerable extent depend upon the evaporation toss. 

The diurnal fluctuation of Tjj(0) fw a dry soil (Figure 42B) was caused primarily 
by soil temperature changes. As a result, the maximum range aT|^( 0) was higher at 
10.0 GHz (20®K) than at 1 ,4 GHz (12°K); this reflects smaller soil temperature fluctuation 
at greater depths (Figure 41). As mentioned before, the 0000 hours and 2400 hours 
brightness temperatures were olmost identical due to small moisture loss (and consequently 
small emissivity change) during that period. Surface configuration had less effect on Tj^(0) 
of dry soil , 

The effect of soil moisture changes on the brighhiess temperature has been 
isolated in Figure 43 which shows total emissivity computed using Equations (22) 
and (23) . Emissivity varied considerably during periods of large evaporation loss 
but remained approximately constant if soil moisture changed little; the latter 
condition existed when either atmospheric evaporative demand was low (nighttime. 

Figure 43A) or soil was dry (Figure 43B). As previously mentioned, the emissivity 
should decrease somewhot during evening hours; such a trend is not shovwi in Figure 
43A. The difference between emissivities at individual frequencies for 0000 hours 
(Figure 43A) was caused by the frequency dependence of the dielectric properties of 
moist soil because soil moisture profiles were identical (Table 14). The same expla- 
nation may be used to account for the time trend of emissivity at 4.0 GHz (Figure 
43A) which increased more during the 24-hour period than those at both 1 ,4 GHz 
and 10.0 GHz; os a consequence of this trend, brightness temperature ot 4.0 GHz 
behaved the same way (Figure 42A) , 

IF all other parameters in Equation (20) remained unchanged, the brightness 
temperature contribution from the layer, AT^. (0) (Equation (21)), would decrease 
with increasing Irequeney because the term exp(" .(0) z.) would decrease. This 

J J 

relationship is illustrated in Figure 44 for the 0-2 cm layer; contributions from 0-1 cm 
and 1-2 cm at 10,0 GHz were added. While for H/0/0 (Figure 44A) only 0,23 of 
^T|^j(0)was contributed from the 0-2 cm depth, the froctlon was 0.75 at 4.0 GHz and 
0,9? at 10.0 GHz (at 10.0 GHz, 0.947 was contributed by the 0-1 cm loyer alone). 

As the soil dried out, contributions from 0-2 cm layers dropped at all three frequencies, 
but the maximum decrease occurred at 4.0 GHz. In contrast, although the dry soil 
^T|jj(0) contribution at 1 .4 GHz decreased to one third of its value for wet soil 
(Figure 44A,B) the absolute change was much smaller (0.23 to 0.08). Consequently, 
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Dimensionless 





Moisl-ure 


Profile 


Depth 

m 

cm 

H/0/0 

1 

l/Q/0 

Tetnperafure at time (hours). In 

00 

06 

12 

18 

24 


00 

06 

12 

18 

24 

0 

20.5 

19.9 

31.1 

27.0 

20.2 


20.0 

19.9 

42.3 

33.5 

20.2 

1 

21.0 

20.1 

30.6 

27.3 

20.8 


20.9 

20.2 

40.9 

34.1 

21.2 

3 

22.0 

20.7 

29.5 

27.9 

21.9 


22.7 

21.0 

38.4 

35.2 

23.1 

5 

23,8 

21.4 

28.7 

28.2 

22.9 


24.4 

21,9 

36.1 

35.7 

24.7 

7 

24.5 

22.1 

27.9 

28.4 

23.9 


25.9 

22.9 

34.2 

35.7 

26.2 

9 

25.2 

22.7 

27.3 

28.5 

24.7 


27.2 

23.9 

32.6 

35.5 

27 .A 

U 

25.9 

23.4 

26.9 

28.4 

25.4 


28.3 

24.8 

31.4 

35.1 

28.5 

13 

26.3 

24.1 

26.5 

28.4 

25.9 


29.4 

25.7 

30.4 

34.6 

29.4 

15 

26.8 

24.6 

26.3 

28.2 

26.4 


30.0 

26.5 

29.7 

34.0 

30.1 

17 

27,1 

25.2 

26.2 

28.1 

26.9 


30.5 

27.3 

29.2 

33.4 

30.6 

19 

.27.4 

25.7 

26.2 

28.0 

27.2 


31.0 

27.9 

29.0 ! 

32.9 : 

31.0 

21 

27.7 

26.1 

26.3 

28.0 

27.5 


31.2 

28.5 

28.8 

32.9 

31.3 

23 

27.9 

26.6 

26.5 

27.9 

27.7 


31.4 

29.0 

28.8 

31.9 

31.4 

25 

28.1 

1 

26.9 

26.7 

27.9 

27.9 


31,5 

29.5 

28.9 

31.6 

31.5 

27 

28.3 

27.3 

26.9 

27.9 

28.1 


31.6 

29.8 

29.1 

31.2 

31.6 

29 

28.4 

27.6 

27.2 

28.0 

28.3 


31.6 

30.2 

29.3 

31.0 

31.6 

31+ 

28.5 

27.9 

27.4 

28.1 

28.4 


31.6 

30.4 

29.5 

30.8 

31.5 
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of time. 
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large moisture changes in the 0"2 cm layer would best be observed at 4,0 GHz (note 
that at 10.0 GHz, the 0-2 cm layer still contributed more than 0.7 of Tj^(O) when 
the soil was dry). In addition to the depth distribution of the contributions, this 
conclusion is supported by the brightness temperature change at 4.0 GHz, For wet 
soil (Figure 42A, 0000 houis), the Tj^(0) equalled that at the lower frequency, while 

in the case of dry soil (Figure 42B), Tjj(0) at 4.0 GHz was similar to that at 10,0 GHz. 
For deeper layers, trends became reversed as the soil dried out. This is illustrated 
in Figure 44 for 4-6 cm and 10-12 cm layers at 1 .4 GHz. It is apparent that the Tbj(0) 
trend for these layers was opposite to that for 0-2 cm, both diurnal ly and for a longer 
period: os the 0-2 cm contribution (and more generally that from shallower layers) 
decreased, that from the deeper layers increased. 

The magnitude of M^j[0) os a function of depth is shown in Figure 45A for 
H/O/O, While the bulk of microwave energy was contributed by the 0-1 cm layer at 
10.0 GHz and 0-4 cm layer at 4.0 GHz, contributions at 1.4 GHz diminished only 
gradually with depth, and some radiation traveled from depth (s) even below 32 cm. 

At all three frequencies, the contributions decreased asymptotically with depth. For a 
dry soil, however, the asymptotical decrease occurred at 10.0 GHz only, while ot the 
lower frequencies the layer of maximum brightness temperature contribution shifted 
downward. This resulted in extending the "effective" depth of penetration at all 
frequencies . 

It is interesting to note that the changes of ^T|^j(0) due to soil drying depended 
on frequency. That is, the ^Tj^j[0) decrease for near^urface layers and consequent 
ATbjO) increase for the deeper layers were confined to a smaller total thickness of the 
soil at the higher frequencies; this resulted in smaller gradients with depth 

at the lower frequencies. For example, between H/O/0 and 1/0/ **^TjjJ(0) decreased 
for the 0-6 cm layers at 1 ,4 GHz, but only for 0-2 cm at 4.0 GHz, and for 0-1 cm 
at 10,0 GHz; beneath these layers, ^Tj^J(0) increased for depths 6-32+ cm (1 ,4 GHz), 
2-14 cm (4,0 GHz), and 1-7 cm (10.0 GHz)^ However, die 6Tj^j(0) increases at 
greater depths were quite small (less than 0.02 at all depths for 1 .4 GHz) and possibly 
negligible. This is shown in Table 15 using Tj^j(0) values for H/0/0 and l/O/O at 
1200 hours. Althot^h the total brightness temperature changed by 65.1^K, only 
contributions from the tc^ 4 cm shifted by more than 2%. Such a small change is 
apparently insufficient to appreciably affect the measured brightness temperature; 
consequently, only the moisture changes closer to the surface can be expected to have 
a marked effect on the measured T|^(0). 
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Figure 45. Fraction of brightness temperature contributed for a wet (A) and a dry (B) soil at three frequencies 
as a function of depth . 



Toble 15. Brightness temperature conft'ibutions for a wet and 
a dry moisture profile as a function of depth at 
1.4 GHz, 1200 hrs. 



Brightness Temperature Contribution 

Depth 

in 


(°K) 

Ty(0) 

(dimensionless) 

cm 

H/0/0 

l/O/O 

H/0/0 

iyo/0 

Difference 

H-L 

(1) 

(2) 

CO 

(4) 

IlgHiil 

HHHII 

0-2 

44.7 

20.2 

0.219 

0.075 

0.144 

2-4 

36.0 

36.6 

0.176 

0.136 

0.040 

4-6 

28.6 

37.6 

0.140 

0.140 

0.000 

6-8 

22.1 

32.2 

0.108 

0.120 

-0.012 

8-10 

17.0 

26.8 

0.083 

0.100 

-0.017 

10-12 

13.1 

22.0 

0.064 

0.082 

-0.018 

12-14 

10.1 

17.9 

0.049 

0.067 

-0.018 

14-16 

7.7 

14.5 

0.038 

0.054 

-0.016 

16-18 

5.9 

11.9 

0.029 

0.044 

-0.015 

18-20 

4.5 

9.7 

0.022 

0.036 

-0.014 

20-22 

3.5 

7.8 

0.017 

0.029 

-0.012 

22->24 

2.7 

6.3 

0.013 

0.024 

-0.011 

24-26 

2.0 

5.1 

0.010 

0.019 

-0.009 

26-28 

1.6 

4.1 

0.008 

0.015 

-0.007 

28-30 

1.2 

3.3 

0.006 

0.013 

-0.007 

30-32 

0.9 

2.8 

0.004 

0.010 

-0.006 

32+ 

2.9 

9.2 

0.014 

0.036 

-0.022 

SUM 

204.5 

268.0 

1.000 

1.000 

0.000 
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While soil tnoisiure contents generally vary with depth, a radiometrically 

measured brightness temperature is an integrated, weighted expression of contributions 

from various depths (Equation (20)). There should exist, then, some single moisture 

value which could be considered representative of that part of the soil profile 

involved in generating the measured microwave radiotion. From Equation (20) it 

would appeor that for a given measured brightness temperature, the driest effective 

moisture estimate is one in whic^ the moisture content of each layer m. is weighted 

lL J 

by the contribution of the j layer T|^j(0) to the total brightness temperature, T|^(0). 

Such a parameter was defined by Equation (24) and designoted m^|^. Colunrn 4 of 

Table 16 contains m^j^ values for H/0/0 and 1/0/0 at the three frequencies for 6-hour 

intervals during the diurnal cycle. As expected, these values (i) changed less than 

and (ii) were higher than near-surface moistures in all cases because of increasing 

moisture contents with depth (Table 13), Since m^|^ values ore affected by the current 

physical temperature of the soil, it is of interest to determine the extent of soil 

temperature influence on the effective soil moisture. This has been done by defining 

tth 


m 


ea 


(Equation (25)) on the basis of emissivity contribU’ion of the ] layer A€j to the 


total emissivity e • The results, presented in column 5 of Table 16, indicate that 
there was little difference between m. values colculated by the two methods; in fact, 
they differed in ail cases by less than 0,40 % of the m^^ value (Table 16, column 
7). The worst case was l/O/O, 4,0 GHz, and 1200 hours with a difference of 0.39%. 

In other words, the various soil temperature profiles had a negligible influence on the 
effective soil moisture computed for various conditions. 

While useful in evaluating the effect of vertical temperature distribution on 

the weighting of cemtributions from individual layers, m and m^j^ do not seem to be 

the moisture contents which are most closely correlated to the computed brightness 

temperatures in the cases studied. Supporting evidence for this statement can be found 

in Figure 46 in which m (Table 16, column 5) is plotted ogainst e (Table 16, 

column 6) for the three frequencies. At 1 ,4 GHz for example, the L/0/0 profile had 

3 3 

emissivity of af^roximately 0.89 although its m^^ was equal to 0.17 cm /cm , Note 
that for a completely dry soil, ^ = 0,91 at this frequency. This means that while 
some water was present within the soil that contributed to the total computed bright- 
ness temperature, the moisture did not substantially affect the value of € . Results 
for the remaining beqjencies were similar, except the curves were shifted toward 
lower moisture contents (Figure 46). In other words, the profiles behaved radiometrically 
almost as if no water were present, although some water always was within the 
contributing depth . 


no 


Table 16. Equivalent moisture contents, emissivity, and the "brightness 
temperature" just below the surface for two profiles os a 
function of frequency and time of day . 



Frequency 

(GHz) 

Time 

Moisture 

(cn?/cm^ 

Emissivity 

{dimen" 

sionless) 


Tb(0/1-) 

(°K) 

Vi 

(°K) 

(hrs.) 



(5) 

. J%)_. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

H/Q/O 

1.4 

00 

0.3000 

0.3000 

0.633 

0.00 

296.8 

294.0 



06 

0.2980 

0.2980 

0.637 

0.00 

295.1 

293.1 



12 

0.2749 

0.2750 

0.678 

- 0.C4 

301.3 

303.6 



18 

0.2474 

0.2474 

0.737 

0.00 

300.2 

300.3 



24 

0.2442 

0.2440 

0.744 

0.08 

295.8 

293.8 


4.0 

00 

0.3000 

0.3000 

0.644 

0.00 

294.4 

294.0 



06 

0.2965 

0.2965 

0.648 

0.00 

293.3 

293.1 



12 

0.2553 

0.2554 

0.694 

-0.04 

303.2 

303.6 



18 

0.2057 

0.2057 

0.760 

0.00 

300.4 

300.3 


- 

24 

0.2002 

0.2001 

0.768 

0.05 

294.2 

293.8 


10.0 

00 

0.3000 

0.3000 

0.669 

0.00 

294.0 

294.0. 



06 

0.2961 

0.2961 

0.671 

0.00 

293.1 

293.1 



12 

0.2498 

6.2498 

0.707 

0.00 

303.6 

303.6 



18 

0.1916 

0.1916 

0.762 

0.00 

300.3 

300.3 



24 

0.1851 

0.1851 

0.768 

0.00 ! 

293.9 

293.8 

i/o/o 

1.4 

00 

0.1693 

0.1687 

0.872 

0.36 ’ 

296.6 

293.9 



06 

0.1691 ! 

0.1686 

0.873 

0.30 

294.1 

293.2 



12 

0.1663 

! 

0.1669 1 

0.875 

“0.36 

303.3 

313.9 



18 

0.1638 

0.1640 

0.878 

"0,12 

304.3 

307.1 



24 

0.1639 

0.1634 

0.878 

0.31 

297.2 

294.2 


4.0 

00 

0.1052 

0.1055 

0.905 

-0.28 

294.5 

293.9 



06 

0.1054 

0.1052 

0.905 

0.19 

292.7 

293.2 



12 

0.1029 

0.1033 

0.906 

-0.39 

309.5 

313.9 



18 

0.1003 

0.1002 

0.908 

0.10 

306.6 

307.1 



24 

OJ0999 

0.0996 

0.908 

0.30 

295.1 

294.2 
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Tflble 16. Equivolenl^ moisture contents, emissivfty, and the “brightness 
temperature” just below the surface for two profiles as a 
function of frequency and time of day. 


Soil 

Profile 

Frequency 

(GHz) 

Time 

Moisture, in (cn?/bm^ | 

Emissivity 

(dimen~ 

sionless) 

(4hg) 

(5) 

_ era) 

3 

Vi 

(°K) 

Oirs.) 


(,(0) 

ea 

1 

2 

3 

4 

5 

6 

7 

8 

9 

l/tV'O 

10.0 

00 

0.0721 

0.0720 

0.915 

0.14 

294.3 

293.9 



06 

0.0720 

0.0720 

0.915 

0.00 

293.1 

293.2 



12 

0.0700 

0.0700 

0.917 

0.00 

312.2 

313.9 



18 

0.0660 

0.0659 

0.922 

0.15 

307.1 

307.1 



24 

0.0657 

0.0656 

0.922 

0.15- 

294.7 

294.2 

. 

t 

j 

j 


r* 


’ 


i 
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To help understand the above results, one con write Equation (20) assuming 
that the reflections at the interfaces between individual layers are negligible. In 
this case, the radiometric temperature which would be measured immediately below 
the air**soii interface ot the normol incidence T^(0, I*) is equal to 

H 

Iffc (0,1-) = S Tp, ( I - exp (- T , (0)Az, )) 

1 i 

. II exp(“ S (0)Az^_,j) . (26) 


For the first three layers, (26) becomes 


+ (Tp2- Tp,)exp(-7,(0)Az,) + ( Tpj - ) 

. exp (“ Y ^(0) Az^ - Y 2 ( 0 ) Az 2 ) 

. Tp 3 exp (- Y ^(0)Az, - Y2(0) Az 2 " Y3(0)Az3 ) . (27) 

From (27) it is envious that if soil temperature is constant with depth and if reflections 
between individual layers are negligible, the radiometric teneperoture which would 
be measured immediately below the soil~air interface is equal to the physical temperature 
of the first layer, . Secondly, if the reflections are negligible but soil temperature 
varies appreciably from layer to layer, T^(0,1^ will also be affected somev^at by 
values and the moishire contents (due to the latter's influence on Y (0)) of the respective 
layers. Finally, if the inrerfdee reflections are not negligible, conhributions from 
individual layers will be further modified by these reflections. 

Values of radiometric ten^roture immediately below the air~soil interface 
Tjj(0, b) were calculated at the three frequencies for H/Q/0 and l/O/O. The results, 
given in Table 16 (column 8) show that especially for H/0/0, T^(0,1') was indeed 
very close to Tp<| (column 9 of Table 16). That is. Internal soil reflections and 
temperature gradients were not iti^ortant in this cose. For the dry soil, larger differences 
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between T|^(0,1") and Tp^ occurred which can be atfributed to the dturnol heat 
wave penetration into the soil. Nevertheless, an inspection of T|^(0, 1~) and 
T ^ results in Table 16 reveals that the maximum difference between the two was 
aCout 4,4*^K which represents only 0,02 of the commonly observed soil temperatures. 

In terms of microwave remote sensing, Tj^(0, 1-) is not the value that is 
measured, however. As expressed by Equation (20), the measured value is Tj^(O) 
which is related to Tjj(0,l“) os follows; 


Tj^(0) = (0,1-)(1-R), 


(28) 


where R is the power reflection coefficient at the air-soil interface. For a plane 
smooth soil surface, the value of R is a function of dielectric properties of the soil 
which in turn depend prlmorily on soil moisture content (Chapter 1). The above 
results thus imply that under conditions where the brightness temperature model of 
Burke and Paris (1975) is valid and where soil moisture and temperature profiles are 
similar to those used here, it Is only the soil moisture content at and near the soil 
surface which determines the brightness temperature measured by microwave radiometers. 
This also explains why the l/O/O profile exhibited high emissivity although it contained 
an appreciable amount of water within the contributing soil profile section; the reason 
suggested by the obove analysis was the low water content in the first soil layer 
(Table 13). It is notew^thy that in an experimental study of the reloticmship between 
soil moisture and rodiometric temperatures measured by Skylab S194 radiometer 
(frequency 1.4 GHz), Eagleman and Uloby (1974) also found that moisture in 0-2.5 cm 
(4 data sets) or 0-5.0 cm (1 data set) was most closely ccNTelated to the microwave 
measurements. The above conclusion does not contradict the previous discussion of 
contributions from various subsurface layers (Figure 44, 45) because in the case of 
small temperature differences and small reflections between individual loyers, the 
radiative energy transfer proceeds similarly as In a homogeneous medium. 

In extrapolating the brightness temperature model predictions to natural conditions, 
it should be kept in mind that the model neglected surface and volume scattering. 

Since scattering Is a form of reflection, the microwave energy can be scattered only 
if a change in dielectric properties exists between the media on the opposite sides of 
the scattering surface. To get an idea of the magnitude of scattering, consider inter- 








3 3 

faces between three media, namely air, dielectrically oniform soil with 0.05 cm /cm 

3 3* • 

moisture, and dielectricolly uniform soil with 0,25 cm /cm moisture content. Using 
coefficients in Table 12 to estimate k' and computing power reflection coefficient 
between media 1 and 2, R^2' from 



the following results can be obtained at 1 ,4 GHz: 


Medium 1 Medium 2 R ^2 

3 3 

air (assumed dry) soil, 0.05 cm /cm 0.09 

soil, 0,05 cm^cm^soil, 0,25 cm^/cm^ 0,06 

3 3 

air (assumed dry) soil, 0.25 cm /cm 0.31 


3 3 

Thus if soil is relatively dry (0.05 cm /cm ) near the surface ond below 

3 3 

this layer moisture suddenly increases to 0.25 cm /cm , the second reflection 
(0.06) and therefore the higher moisture would have less impact on the microwave 
signal than the first one (0,09). An air-moist soil (0.25 cm /cm ) interface would 
have a higher effect, however (R^2 ~ 0.31). These considerations seem to suggest 
that subsurface scattering (and therefore the effect of subsurface moisture on microwave 
signal) will be higher for a soil consisting of more or less homogeneous clods separated 
by air voids; these clods should be larger than approximately one tenth of the wove* 
length measured in order for scattering to take place. On the other hand, in a 
relatively homogeneous soil (in terrm of air*soil material distribution), even sudden 
large moisture changes do not seem to have a major effect on the microwave radiation 
(R|2 ~ 0.06) and ther«$fore their impact on the measured signal will be smoller. 
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3 3 

It should also be noted Hiat the assumed soil moisture chonge frwn 0,05 cm /cm 
3 3 

to 0.25 cm /cm is unrealistically large compared to coses which might be en- 

countered under field conditions. For example, the maximum moisture change 

with depth reported by Jackson (1973) for the July, 1970 experiment (average 

daily potential evaporation of 1.077 cm) wos from 0.079 (depth 0“1 cm) to 

0,168 cmVcm^ (1“2 cm), or a change of approximately 0.09 crt^/cm^ compared 
3 3 

to 0,20 cm /cm*^ used here. Furthermore, the radiation transfer in the soil proceeds 
continuously and not in discrete steps, and soil moisture also changes in a more or 
less continuous foshion with increasing depth; therefore, the chance of encountering 
a sudden moisture change oppears relatively small. 

If only the surface and near~surface moisture contents ore of major importance 
in determining' the microwave brightness temperature, then it follows that the soil 
surfoce temperature may be used as the representative soil temperature without 
generating appreciable errors. The passive microwave sensois could thus be supple- 
mented with infrared radiometers to yield soil emissivities (computed from Equation 
(2)), This sensor combinati<x> was proposed by Poe et al. (1971), Consequently, 
predictions of soil temperature would not be necessary provided that physical soil 
temperature values can be isoloted from the infrared measurements for areas correspond- 
ing to the microwave resolution cells. 

Brightness temperature calculations discussed so far in section 4.1 were made 
with the assumption that dielectric properties of a moist soil depend on water content 
only. This is because the temperature dependence of the soil dielectric constant 
has not been studied experimentally, most laboratory measurements having been made 
at approximately 20°C. On the other hand, the dependence of the dielectric 
constant of woter on temperature has been well established (Peris, 1969), To deter- 
mine the order of magnitude for temperature effect on dielectric properties of a moist 
soil, the following procedure was used: 

(i) Calculation of the ratio between dielectric constant of water at 
5°C and 35°C and that at 20°C. This was done separately for k' 
and k", at each frequency, and for salinities 0 *>6o, 10 ^oo, and 
20 Voo. To obtain the ratios for k* and k” as a function of 
frequency and temperature, values for the different salinities were 
averaged. 
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(ii) Using ihese ratios and dielectric constant values of a moist 
loam (Cihlar and Ulaby, 1974), k' and k*' for a given 
temperature (5°C or 50®C) were calculated by multiplying 
the measured values by the above ratio and the volumetric 
water content. In other words, the assumption was made that 
the change of the dielectric constant of a moist soil with 
temperature is solely due to the temperature dependence of the 
dielectric constant of the water portion. Except for possible 
ionic complex interactions, such an assumption oppears 
justified since the dielectric properties of dry soil are 
practically constant with temperature (Cihlar and Ulaby, 1974). 

(iii) Power reflection coefficient and emissivity were computed for 

soil profiles in which moisture was constant with depth (profiles 
3 3 3 3 

between 0.00 cm /cm and 0.35 cm /cm ) and temperature was 
also constant (either 5°C or 50°C), 

Figure 47 shows the results at three frequencies; all 5*^C data are joined by 

curves. Apart from the frequency dependence, these data also show a systematic 

effect of temperature on power reflection coefficient and emissivity: at all three 

frequencies, the power reflection coefficient decreased with increasing frequency, 

particularly at higher moisture content. The maximum difference (0.007) corresponds 

3 3 

to a moisture content difference of approximately 0.010 cm /cm . Considering the 
relatively extreme soil temperatures which this moisture differential represents, it 
would appear that the temperature effect on the soil dielectric properties could be 
neglectfcd in an algorithm designed to extract soil moisture information from microwave 
measurements without introducing appreciable errors. 

As stoted in section 4.1, calculations of brightness temperatures T^(0) (and 
therefore effective emissivities e (0)) were made for incoherent radiation. Active 
sensors transmit a coherent beam, however, and the question thus arises whether the 
same relation between power reflection coefficient R and emissivity g- , e = 1 “ R, 
is valid for both coherent and incoherent radiation. To answer this question, power 
reflecticxi coefficients were calculated for 11 soil moisture profiles in two ways: 

(i) R = 1 - e 

where e = e (0) from the brightness temperature model (Equation (23)), 
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Power Reflection Coefficient, Dimensionless 














(li) R = 

Values of p were calcuiat'ed* by a method outlined by Casey (1973) 
in which the second-order linear differential equations describing the 
electromagnetic field behavior in the vertical direction within the soil 
are solved in terms of Hill's functions. 

A comparison of results obtained by the two methods is given in Table 17, 

The two power reflection coefficients are almost identical at higher frequencies and 
high moisture contents, but the agreement is also very close at 1 ,4 GHz if the moisture 
content is high. Greatest discrepancies occurred at 1 .4 GHz when moisture content 
was low. The discrepancy between power reflection coefficients computed for coherent 
and incoherent radiation thus appeared to increase with increasing depth of penetration 
or vertical nonuniformity of soil moisture contents (see Table 13} or both. This agree- 
ment also suggests thot under conditions when the assumptions of the plane Interfaci^, 
multi-layer, nonscattering soil material ore valid, the return of coherent active mi .*o- 
wave signal will be determined by the surface and near surface moisture content, 
similarly as was the case with brightness temperature , In addition, the discussion of 
diurnal variation of emissivity may be directly extrapolated to power reflection coefficient. 

4.2 SOIL MOISTURE CONTENT DETERMINATION OVER LARGE AREAS 

The following discussion is based on two premises. First, soil water coritent 

should be specified in volumetric units in order to provide usable information, be it 

3 3 

equivalent depth (in cm), amount of water per unit soil volume (in cm /cm ), or 
some other parameter. Secondly, the soil moisture results should be location- 
specific to be of maximum benefit to the user,. While the first premise rests on the 
fact that water is involved in various processes and environments as quantities, 
the second assumption appears Justified by the grov. ng environmental concern about 
conditions at specific locations (Holton and Lopez, 1970), as well as economics of 
land use of particulor areas (such as fields). 

Results of the previous section and consideration of the effect of moisture on 
the electromagnetic radiation propagation suggest that microwave remote sensors may 

*Calculations were made by Dr.K.F. Casey, Department of Electrical Engineering, 

Kansas State University, Manhattan, Kansas. 
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Frequency (GHz) 


Time 

1, 

.4 

4.0 


Bi 

I* 

C 

-L - 

06 

0.364 

0.363 

0.349 

0.352 

12 

0.306 

0.322 

0.299 

0.306 

18 

0.250 

0.262 

0.230 

0.240 

24 

0.246 

0.256 

0.223 

0.232 

18 

0.251 

0.264 

0.232 

0.241 

18 

0.256 

0.269 

0.239 

0.249 

18 

0.255 

0.267 

0.237 

0.246 

06 

0.073 

0.118 

0.053 j 

0.082 

12 1 

0.072 

0.113 

0.050 

0.078 

18 

0.072 

0.109 

0.047 

0.074 

24 

0*072 

0.109 

0.048 

0.074 


*Nofe C = coherent case I = incoherent cose 



0.327 

0.290 

0.233 

0.226 

0.234 

0.241 

0.239 

0.047 

N.A. 

N.A. 

















provide soil moisi’ure Informai’ioh abouf the soil surface and near subsurface layers 
unless the soil material appears geometrically or dielectrically rough relative to 
the wavelength measured. Since various types of land use and management are 
concerned with deeper layers as well(e.g. profile thickness of 100 to 200 cm ), it 
follows that a microwave remote sensing technique of soil moisture estimation would 
hove to become a part of a larger system involving other inputs. Following is a 
discussion of a possible algorithm of which the objective is providing slte'^pecific 
profile soil moisture content information during frost-free periods over areas with 
agriculture os the main land use category. 

The total area of Interest is divided into a number of small cells i of equal 
size, thus forming an equal area grid system (Tomlinson, 1972). The maximum cell 
size must be smaller than the smallest area element the soil water content of which 
is needed and larger than the minimum resolution cell size of the microwave sensor. 
The algorithm is based on the soil water balance equation (Equation (4), (7)} which, in 
the absence of irrigation and after 's replaced by -(R4D), may be written as 


S(t+At) = S(t) + P(At) - R(At) - AE(At) - D (At), (29) 


where R - surface runoff. To implement Equation (29) in a site-specific algorithm, 
the following information is needed for each area element i and time increment At: 
precipitation, wnoff, actual evaporation or actual evapotranspiraticm, and drainage 
below the terminal depth of the profile. If the remote sensing mission takes place 
every other doy, At = 2, and the remaining necessary parameters may be obtained 
for the period of At as follows. 
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PrecipUart<»i 

Precipi^aHon amounts P ( A t)and duration PD can be acquired from 
meteorological statioru as values representative of the cells in which they are 
located. Point-rainfall extrapolation methods (Hutchinson and Wal ley, 1972; 
Chidley and Keys, 197Q; Unwinn, 1969) can then be used to obtain areal rainfall . 

As mentioned in Chapter 2, however, precipitation patterns vary considerably even 
at short distances, and areal rainfall may be far from representative for specific 
locaticms. To solve this problem, microwave sensors can be used for detennining the 
actual precipitation patterns over the entire areo. Thus, the actual precipitation 
for the cell i con be considered a) equal to the areal rainlbll computed above if 
the microwave remote sensing data show that cell i received precipitation, b) equal 
to ze:o otherwise . 

Runoff 

Runoff can be computed as the difference between precipitation Intensity and 
infiltration capacity. If the infiltration capacity is calculated by Holtan's (1970) 
formula, runoff R( At) for cell i and precipitation (P( At)), is equal to 


R (At) = P (At) - PD ( f + a ), (30) 

c o 

where 

a = infiltration capacity in inches/hour/(inch) ^ of available storage; 

= available storage in the surface layer (*'A“ horizon in c^ricultural 

soils) in inches of equivalent water Icyei; 

f = constant rate of infiltration after prolonged wetting in inches/hour, 
c 

Values of a were given by Holtan et ol . (1974) for various crop types and thc^e of f^ 
by Musgrave (1955) for different soil hydrologic classes. Considering possible multiple 
precipitation events as one will not result in large errors because of At (Richards 
and Strahl, 1969). The runoff R(At) can be either routed to the channel (Holtan 
et ol ., 1974) or simply ignored; the former strategy would help to account for that 
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pa^^ of R( which enters the soil prior to its reaching the channel. At this point, 
the actual amount that entered the soil, At), can be calculated for each cell I 
as 

P^ff (At) = P (At) - R(At) . (31) 


Actual EvgporatiiMi or Evapotranspiration 

Various methods can be used to compute actual evaporation (section 2. •). 

Hoitan et al. (1974) described a technique for evapotranspiration estimation which 
is based on pan^evaporation data, growth index, and soil moisture status. 

Drainage 

The drainage component D( At) can be estimated with the assumption that only 

a certain maximum amount of water, S^, con be retained for extended periods of time 

(section 2.3). Results in section 2.3 olso suggest thot D( At) may be considered 

negligible if the soil storage capacity S is not exceeded. Consequently, droinoge 

P 

can be determined as 

D(At) = 0 if Pgff(At)^ Sp - S(t) 

**eff (At) + S (t) “ Sp otherwise. (32) 

A more refined approach to computing D( At) wwld consist of combining Equations (32) 
and (30). 

The values of P, R, AE, and D must be determined for each cell i and period At. 
Then, given the initial soil moisture storage S(t), the amount of water stored at time 
(t + At) can be computed from Equation (29). Figure 48 is a flowchart of the sequence 
of computations for one At period. 

The procedure described above would yield one moisture storage value for the 
entire soil profile. As pointed out in Chapter 2, however, it is preferrable to subdivide 
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Figure 48. Flowchart of an algorithm for large area monitoring of soil 
water content. 







fhe profile info zones and monifor moisfure changes for each zone separately. For 
a bare soil, the two zones model used in Chopter 2 can be used with the upper zone 
divided into layers as thin os 1 cm if such resolution is warranted; the lower zone 
would comprise the remainder of the profile (section 2.4). In the case of vegetated 
areas, the loyered model of Baier et al. (1972) would be more appropriate because 
it allows for changing root structure of the plants. In the multi**tayer case, drainage 
from one loyer can be regarded as input into the layer immediately below. 

To inclement the above algorithm for monitoring soil water content, the 
following data would be needed; 

a) Time ** invariant parameters: 

~ Hydrological soil group to determine f (Equation (30)). 

- Storage capacities of individual zanes within the 

soil profile; these could be cbtoined from soil survey reports. 

b) Time - varying parameters: 

- Land use (crop type), to be specified for longer periods (e.g., 
a season) . 

" Precipitotion amounts, durations, and meteorological variables 
needed for AE ( t) calculotions. 

The above described algorithm would allow combining advantc^es of various 
methods of soil moisture estimation useful for large areas (section 1.1) while avoiding 
their drawbacks. Furtheitnore, it would employ data readily available in many areas. 
It should be noted that vari<»js hydrological models (Knapp, 1973; Staff, Hydrological 
Research Laboratory, 1972; Holton et al., 1974) approach soil moisture accounting in 
a way similar to that used here. The difference is that soil moisture was the only 
output parameter required which allowed simultaneous simplification and exponsion 
compared to the other models. 
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CHAPTERS. 


MICROWAVE RESPONSE TO SOIL MOISTURE: A REVIEW 
OF EXPERIMENTAL MEASUREMENTS 

As stafed in Chapi’er 1, I’he concept’ of soil moisture estimation from remote 
microwave sensor measurements is based on the microwave signal response to the 
presence of water in the soil . This phenomenon has been demonstrated in Chapter 1 
(Figure 2) and Chapter 4 by applying simplified theoretical solutions to vertical soil 
temperature and/or moisture profiles which were either constant with depth or 
computed from empirical models . In order to investigate the role of moisture and 
temperature regimes in microwave remote sensing of soil moisture/ only the intrinsic 
relationships between microwave reflection (and emission) and soil moisture content 
were considered. Under natural conditions, however, other target parameters also 
effect the scattering and emission properties of the surface; namely, roughness and 
vegetation cover. 

The difficulty involved in relating some of the statistical parameters used in 
theoretical scattering and emission models to easily specifiable target parameters has 
led investigators to construct empirical models based on expenmental data. One of 
the major objectives of microwave remote sensing of soil moisture programs is to estab~ 
lish sensor parameters (frequency, polarization and ar^le of incidence rcmge) for which 
ti^ effects of target characteristics other than soil moistuie are at a minimum. Hie 
dsjective of this chapter is to present a summary of active and passive microwave 
measurement progroms in remote sensing of soil moisture. The presentation is divided 
into sections according to platform (ground, airborne, spaceborne) and each section 
is divided into active and passive observations. 

5.1 GROUND BASED PLATFORMS 

5.1.1 Active Microwove Observotions 

Backscatter data acquired by a truck mounted Active Microwave Spectrometer 
over the 4-8 GHz frequency region were reported by Uloby (1974) for a bare field. 
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The soil type was Pawnee clay loam and the surface was slightly rough (2.5 cm rms 
he^ht). With soil moisture represented by the average gravimetric moisture content 
in the top 5 cm of the soil, m^ (0,5)*, the backscattering coefficient showed a 
linear increase with soil moisture after the latter quantity had exceeded 15%; no change 
in o® was observed for moisture contents below 15%, 

In an attempt to provide a more representative value for soil moisture, Ulaby 
et al . (1974) computed an "effective" moisture content, m^(0, S ) in grams per cm , 
where m^(0, $ ) was defined as the average moisture content within the skin depth . 

Plots of os a function of m^(0,5) and m^(0, 5 ) are shovm in Figures 49 and 50, 
respectively. With this new approach for defining moisture content, o° indicates o 
linear response over the entire range of m^(0,S). The radar signal change with soil 
moisture was highest at low incidence angles (Figure 50) and the highest frequency 
(7.1 GHz) and decreased as the incidence angle increased from 0^ to 70^ and as 
frequency decreased from 7.1 GHz to 4.7 GHz . These trends were quantified in terms 
of sensitivity, defined as the change in scattering coefficient resulting from a change in 
m (0, j) of 0.01 g/cm (Figure 51). The various frequency/polarization combinations 
exhibited an approximately equal sensitivity at 20^ incidence angle for the Pawnee 
cloy loam data. 

Because of the range of roughness choracteristics that exist under field 
conditions, the most reliable method for determining the effect of roughness on the 
sensitivity of microwave sensors to soil moisture variations is by analysis of experimental 
data. During the summer of 1974, o° data were acquired for each of three bare 
fields having identical soil type but different scales of surface roughness ( Batlivola 
and Ulaby, 1975). The three fields, which for identification purposes will be refen^ed 
to as smooth, medium rough and rough, hod rms surface heights of 0.88 cm, 2.6 cm 
and 4.3 cm respectively. Simultaneous with the radar measurements, passive radiometric 
measurements were also acquired by Texas A&M University. Analysis of rtie passive data, 
however, has not yet been completed. The active microwave measurements were 
conducted at 8 frequencies between 2 and 8 GHz, for all linear polarization configurations 
at angles of incidence between nadir and 40^ ( in 10^ increments). 


* In Chapter 5, the parameter used to describe soil woter content is m^^,c), where a 
is w (moisture in % by weight) or v (moisture in g/cm^), and b and c represent upper 
and lower depths of tne layer represented (in cm) . 
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Figure 49. Scattering coefficient as a function of moisture ■ J 

content in the top 5 cm. Frequency is 4. 7 GHz. I ^ 


From Uloby, Cihlar and Moore (1^4)_ 
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Figure 50. Scattering coefficient as a function of effective 
moisture content. Frequency is 4. 7 GHz. 


From Uioby, Cihlor and Moore (1974), 









Figure 51. MoiStUre 




After acquiring radar data with the fields dry, the three fields were sprinkled 
with water until saturation. Each of the fields was then periodically monitored with 
active microwave measurements and ground truth equipment as the soil dried up. 

Figure 52 Is a time plot of moisture content and cr*^ at nadir, 10° and 20° for 
VV polarization at 2.75 GHz. The data are for the medium rough field. Based on 
moisture samples acquired at 8 sampling sites on the field, the mean and standard 
deviotion of the moisture content were calculated; the moisture Interval shown in 
Figure 52 represents the mean + standard deviotion. The large uncertainty in soil 
moisture content, represented by the width of this interval,! I lustrates the difficulty in 
correlating o“° with moisture content. The uncertainty in m^(0, 6 ) is attributed in 
part to the difficulty in sampling thin loyers of soil and in part to the spatial variatitm 
of moisture content across a field. 

Comparison of cr° to m^(0, S ) shows good correlation at all angles. Similar 
results ai% obtained at other frequencies and polarizations for each of the three fields 
individually. To demonstrate the effect of surface roughness. Figure 53 shows the angular 
response of for the three fields at approximately the same moisture content. Vi^ereos 
the smooth field at 2.75 GHz shows a drop of about 23 dB between nadir and 10°, ifie 
rough field, on Hie other extreme, shows less than 1 dB change. The interesting 
observation in Figure 53 is that at each of the three frequencies shown, a narrow range 
in angle of incidence exists for which the effect of roughness is at a minimutp; at 5.25 GHz 
10° appears as a crossover point and at 7.25 GHz,20° is the crossover point. Appli- 
cation of this phenomenon is shown in Figure 54 where o'® is plotted as a function of 
m^(0,<5 ) for all three fields combined. The frequency is 4.75 GHz, the angle of 
incidence is 10° and the polarization is W. Even with the uncertainty associated with 
the values of m^ (0,6), the correlation coefficient is 0.69. 

Vegetation cover hinders the microwave remote sensing of soil moisture. In 
oddition to attenuating the backscattered or emitted soil radiation, the vegetation 
canopy adds its own contribution to the radiation measured by the radar or radiometer. 

The degree to which vegetation affects the capabilities of microwave sensors in soil 
moisture content determination is a function of plant geometry and the dielectric properties 
of the canopy volume, both of which are related to plant type, stage of growth, stand- 
ing biomass, and other soil and plant variables. Given the complexity of inferoctions 
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Figure 54, Sca^^ering coefficieni’ response as a funcHon of sol! moisfrure 

for three surface roughness profiles combined. From Uloby, et of. (1975), 


among these factors and the lack of understanding of how they relate to one another 
and to the microwave response of vegetation , experimental measurements of the 
soil/canopy system under a variety of conditions are the most efficient method for 
determining the effects of plant canopy on microwave measurement of soil moisture. 

So far, the attention has been focused primarily on agricultural crops. 

Figure 55 (Ulaby, 1975) illustrates the effect of frequency on the rador 
sensitivity to soil moisture differences for a field of mature milo. At nadir, a large 
separation is observed between the frequency responses of dry (8.2%) and wet (26.8%) 
soils. This separation, however, decreases by about 7 dB between 4 and 8 GHz. In 
contrast, a much smaller difference is observed at 30^, particularly at the higher 
frequencies . 

Backscatter data from 43 fields of mature com, milo, soybeans and alfalfa are 
shown in Figure 56 where a® data at 10*^ and 30° are plotted as a function of gravimetric 
moisture content in the top 5 cm (Ulaby, 1975 ). The frequency and polarization ore 
4.7 GHz and HH, respectively. The calculated sensitivity (slope) at 10° is 0.3 dB/1% 
moisture. Similar calculations were performed at other frequencies and armies of 
incidence and are summarized in Figure 57. 

An optimum active sensor for microwave remote sensing of soil moisture should 
operate at the frequency, incidence angle, and polarization for which (i) the cmrrelatioh 
between a° and soil moisture is very high, thus implying low interference by other 
variables, and (ii) the <7° change with changing moisture content is also very high, 
thereby enabling accurate moisture estimates. Figure 57b indicates that practically 
at all incidence angles, the conrelation between and moisture was higher at the 
low frequency, and the difference between polorizatiorts was small. Furthermore, the 
sensitivity S^, was highest at 0° and 10° but decreased quite rapidly for higher 
incidence angles, particularly at 7.1 GHz (Figure 57a). Consequently, Ulaby*s 
(1975 ) study suggested that the low frequency and an incidence angle around 10° 
would exhibit highest moisture sensitivity for the ranges of sensor parameters and 
experimental conditions investigated. 

5.1.2 Possive Microwove Observations 

The passive microwave response to soil moisture has been demonstrated by several 
investigators (Poe, 1971; Edgerton, 1968; Edgerton et al., 1968; Blinnand Quade, 

1973; Lee, 1974; Newton etal., 1974). An example is shown in Figure 58 (Blinn 
and Quade, 1973) where the measured emissivity is plotted as a function of moisture 



Crop Type 


Polarizatio 


Incidence Angle 


Average 

Moisture 

Desiqnati 


0® 

0® 

30” 

30“ 

8.2* 

26.8% 

8.2* 

26.8% 








C C 
0> P 


o 0 


i -5 

to 


43 47 5.1 5.5 5.9 6.3 6.7 7.1 7.5 7.8 
Frequency in GHz 

Figure 55. 4-8 GHz SPECTRAL RESPONSE OF MILO AT o = 0® 

AND 30” FOR LOW AND HIGH SOIL MOISTURE CONDITIONS. 

From Uloby (1975). 











Scattering Coefficient (f In dB 


Crop Height Freq. Incidence 
Type (m) GHz Pol. Angle 


Corn 2.4 4.7 HH • 10® 



Percent Moisture Content by Weight 


Figure 56 . Scal-tering coefficienr response cb a funcHon of soil inoisfure 

for enfire dofa sef regardless of crop fype. From Ufaby (1975). 


Frequency : 4. 7 GHz 
Polarization : HK 



I ncidence Angle in Degrees 


j 1 



Figure lOa. HH Polarization 


Figure Ifflj. W Polarization 


Frequency : 4.7 GHz 
Polarization : Cross 


Crop Type 

Corn 

Milo 

Soybeans 

Altaic 

Crop Height 
in Meters 

2.4 

1.0 

1.0 

as 

Average Soil 
Moisture 

10.6% 

&2% 

8.1% 

10.5% 

Designation 





b 

• 




10 20 30 40 50 60 70 

Incidence Angle in Degrees 


Figure 10c. Cross Polarization 


Figure 57. SCATTERING COEFFICIENT RESPONSE AT 4. 7 GHz FOR 
FOUR CROP TYPES. SOIL MOI STURE CONDITION I S 
LOW. 


From Uloby (1975) 



(o) 


> 

c7) 

v% 


1.0 

A \ ^ i ^ \ ' 

(c) 1.0 

II 1 1 1 I 




» 1 1 

0.9 

_ O _ 

0.9 



a 


8 • 

0,8 

a 

^ 0.8 
cn 

m ^ 

D # 


• 

t/i 

D 

0.7 

■ CAVE ROCK a ^ 


GODECKE ® 4 


■ LOAMY SAND 


CLAY loam. a p. 

0.6 

- •. 10.6W 

0.6 

• 10.6W 


o 10.69H •- 


□ 10.69H 

0.5 

■ 

0-5 

. 1 . 


1.0 

0.9 

£ 

3: 0.8 

in 

j/> 

S 0.7 
0.6 
0.5l 


0 5 10 IS 20 

MOISTUfff CONTENT, % 

1 ' Jn ' I ‘ ' 

-O 00 o fl 


: 


0 


8 


g ^ 


CAVE ROCK 
LOAMY SAND • 
o 31.4V • 1.42V ■ 8-! 
I- a 31.4H ■ 1.42H 
■ ‘ I < I 

5 10 1., cw 

MOISTURE CONTENT. % 


5 10 15 20 25 

MOISTURE C(M0tT, % 


(d) 


«/• 

m 



■ 1 » ■ 

T — 1 1 -T - I 


l.O 


8 Q> 



0 



0.9 

r 1 

• t ® 


0.8 

- 

1 

8 - 

0.7 

- 


SB 

GODECKE 



CLAY LOAM _ 

• 

0.6 

- o 31.4V 

• 1.42V ■ 

i - 


• a 31.4H 

■ 1.42H 

1 

0-5 

■ 1 ■ 

1 1 1 1 1 

1 


5 10 15 20 25 

MOISTURE CONTENT, » 


Figure 58. Moisture effects for loamy sand and cloy loam. 
From Blinn and Quade (1973). 
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Figure 59. Moisture and roughness effects for sandy loom. 
From Blinn and Quade (1973). 
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content' at 1 .42, 10.69 and 31 .4 GHz. The angle of incidence is 25° from nadir. 

Blinn and Quade (1973) reported that the emissivity remainded constant for moisture 
contents below a "breakpoint" ranging between 1.5% and 15%. The breakpoint 
increased with frequency and with the soil cloy content. The moisture content is 
presumably the average moisture content over a 15 cm layer. The existence of such a 
breakpoint is similar to the radar response shown earlier in Figure 49. Moreover, 
the explanation for the change in the value of the breakpoint between 1 .5% and 15% 
as a function of frequency and clay content is in agreement with the change in linearity 
observed for the active microwave data when moisture content was expressed in terms 
of an effective value over the skin depth. The effective layer contributing to the 
emission is closer to 15 cm at 1 .42 GHz than at 31 .4 GHz . 

In addition to measuring the emissivity as a function of moisture content, Blinn 
and Quade (1973) investigoted the effect of surface roughness on the emissivity. 

Regular roughness patterns were prepared by pulling a triangular toothed template over 
the surface of the soil . Tooth spacings of 2.5, 5, 7.5 and 10 cm were used. The grooves 
Were oriented first in the direction of horizontal polorization and then in the direction 
of vertical polarization. Their results, shown in Figure 59/ indicate that increasing 
roughness (furrow spacing) can have a substantial effect on the emissivity, particularly 
at the higher frequencies. 

The effects of surface roughness and vegetation cover were investigated by 
Texas A&M University (Lee, 1974; Newton etai., 1974). Radiometric measurements were 
ocquired at L"band (1 .42 GHz) and X-band (10.69 GHz) with horizontal and vertical po" 
larizotions. Three fields with surface roughnesses comparable to those shown in Figure 53 
were used. With regard to surface roughness, Lee's (1974) results confirm the findings 
of Blinn and Quade (1973) as indicated by Figure 60. The effect of vegetation is illus" 
trated in Figures 61a-c where brightness temperature values of oatgrass can be compared, 
to the response observed with the fields bare. All the data shown in Figure 61 is at 30° 
origle of incidence, L^band. The vegetation appears to produce lower brightness 
temperoture volues than those obtained for the bare soil, particularly for the medium 
rough and rough fields. Similar results were also evident in data at other angles of 
incidence and for different definitions of moisture content (0*18 cm depth and equivalent 
soil moisture). At X*band, the vegetation appeared to completely mask the radiometric 
sensitivity to soil moisturt . 
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5.2 AIRBORNE PLATFORMS 

5.2.1 Active Microwave Observations 

Using airborne scatterometer measurements, Dickey et al , (1974) analyzed 
radar return from dry and irrigated portions of the same fields. The crops grown 
included com, alfalfa, grain sorghum, sugar beets and hay. They found that at 
13.3 GHz and for incidence angles below approximately 45°, values measured 
from irrigated sections were 5 to 7 dB higher than values from the dry sections; 
an example for a corn field (48 cm tall, coverage 40%) is shown in Figure 62. The 
difference in angular responses varied with crop type, however. For example, in 
one alfalfa field, c® for the wet section was only 2 dB higher at angles obove 20°. 

The response to moisture was also observed when all dry section (4239 resolution cells) 
and Wet section measurements (85 cells) were combined (Figure 63). Dickey et al. 
(1974) concluded that similar effects occurred at 0.4 GHz but they were obscured 
because of the large resolution element size. Ulaby et al . (1975) utilized some data 
from the same mission to illustrate the relationship between 15) and at 

13.3 GHz for seven fields with various covers (Figure 64). In spite of the variety of 
crops, cP increased nearly linearly with increasing moisture at 10°. The response 
at 20° was weaker except for one moist field; this anomaly could be related to row 
configuration . At 50°, a° did not change with moisture content in a monotonous 
fashion; this is consistent with the data discussed in section 5.1.1. 

5.2.2 Passive Microwave Observations 

Several airborne experiments have been conducted to determine the effects of 
soil moisture on the microwave emission (Poe and Edgerton, 1971; Jean, 1971; 
Schmugge et al ., 1974). In most cases, an aggregate soil moisture sample from the 
top 15 cm layer was used for defining moisture content. 

Brightness temperature data acquired at 19.35 GHz over an agricultural test 
site in the vicinity of Phoenix, Arizona, was reported by Schmugge et al. (1974) as 
a function of moisture content in the top 15 cm. A breakpoint of about 15% moisture 
followed by a linear decrease of brightness temperature with moisture was observed, in 
agreement with Blinn and Quade's (1973) findings (Figure 59), Data acquired at 
1 .4 GHz for the same fields, however, showed a similar linear decrease but no clearly 
discernible breakpoint was apparent. Schmugge et al. (1974) attributed the breakpoint 
effect in part to the behavior of the dielectric properties of soil as a function of 


142 


-J 


m&lshire contenf', and in part to the Inappropriate representation of moisture content 
by the 15 cm layer. Whereas the average moisture in a 15 cm ioyer may not vary 
much between different fields, the moistures in the top surface (oyer may be 
considerably different. Differences in skin depth between the 19,35 GHz ( X = i.55 
cm) and 1 .4 GHz ( X = 21 cm) frequencies moy partially explain the absence of a 
breakpoint in the 1 .4 GHz brightness temperature response to moisture . 

In a more recent analysis of the data by Schmugge (in; Uloby et al ., 1975), 
the effect of soil type was accounted for by plotting the data for heavy soils (cloy 
loam) and light soils (sandy loam and loam) separately as shown in Figure 65 . Also, 
the moisture content is that of the top 1 cm. The differences in response to soil moishire 
(slope) was then eliminated by plotting the combined data as a iimction of field 
capacity (Figure 66). The result is on improvement in the correlation coefficient over 
either of the plots shown in Figure 65 . Similar results were olso obtained for the 21 cm 
data as shown in Figure 67 . 

5.3 SATELLITE PLATFORMS 


Microwave observations of the earth hove been conducted by both Nimbus 5 
and Skylab sensors. Aboard Nimbus 5 was flown the 19.35 GHz Electrically Scanning 
Microwave Radiometer (ESMR) and Skylab carried a 13.9 GHz Radiometer-Scatterometer 
(S“193) and a 1.4 GHz nadir-looking radiometer (S“194). 

During Skylab's June 5, 1973 pass over a test site in Texas, soil samples were 
collected by ground crews and later anolyzed for their moisture contents (Eagleman, 
1974), The test site was about 300 km in length, which for the S“194 footprint of 
approximately 120 km in diameter, corresponds to only 3 independent data points. 
Successive data points contain about 94% overlap of the area seen on the ground. 
Brightness temperature plots against soil moisture (Eagleman, 1974) produce a correlation 
coefficient of **0.97 with oil the data included. Since inherent in the data is at 
least 94% correlation due to the overlaps, it is not possible to determine the correlation 
with moisture content alone. Hence, S**194 data will not be reproduced in this 
report. 

At 13.9 GHz, S-193 footprints were about 10 times smaller than S-194, 
thereby making it pcssible to investigate the sensitivities of the radiometer and the 
radar to soil moisture content. Figures 68 and 69 ore plots of the emissivity and the 
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Ulaby, etal. (1975). 
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Figure 64. Scattering coefficient as a function of soil moisture content 
at 13.3 GHz, W polarization. Garden City, Kansas. 

From Ulaby et o| . (1975), 
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Figure 68. Jhe emissivity as a function of soil moisture content 

determined from Skylab SI 93 data (2,1 cm wavelength). 

From Ec^teman and Uloby (1974). 
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Figure 69. The SI 93 backseatter coefficient os a function of soil moisture 

for the June 5, 1973 EREP pass. From Eogleman and Ulaby (1974) 




scalt’ering eoefficienl* derived from S~193 meosuremenis over the Texas site 
(Eogiemon and Uioby, 1974). The moisture content value of each of the points 
shewn in the plots is an average of 2 to 3 soil samples gathered within the 
c^proximately 100 km footprint. Hence, some of the points may not necessarily 
be good representatives of the moisture content ofthe entire footprint. Under these 
circumstances, the correlations shown in Figures 68 and 69 ore indeed very encouroig‘‘ 
ing. 
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CHAPTER6 . 


SUMMARY AND CONCLUSIONS 

This study was approached on the premise that an understanding of problems 
Inherent in the mlcrowove approach to soil water content determlnotlon as well cb 
the development of an operationally efficient algorithm for an extraction of moisture 
Information from the microwave data will not be completely successful unless 
consideration Is also given to the major parameters of the soil Involved, namely its 
moisture and temperature regimes . This necessitated a rather brood-based research 
encompassing problems traditionally approached in several scientific disciplines. 

To render the w<»’k manageable, mathematical models were employed wherever 
possible. Furthermore, the scqse was limited to several specific problems (section 
1 .2) and a set of environmental conditions, the major of which were bare soil in an 
agricultural area, cloudless sky, high radiation regime, deep ground water level, and 
absence of scattering at the microwave frequencies of interest. However, care was 
taken to ensure that the mathematical models reproduced reality reasonably well so 
that data generated by the models could be considered representative of the actual 
conditions. Soil moisture and temperature profiles generated by the models for 
various conditions were analy7ed at three frequencies (1 .4 GHz, 4.0 GHz, 10.0 
GHz) using a brightness temperature model and, for a portion of the data,using a 
coherent radiation model which computed an effective power reflection coefficient. 

It should be stressed that the choice of frequencies and incidence angle (0*^) and 
the emphasis on passive mode represent on attempt to accurately reflect the effect of 
moisture and temperature regimes on the microwave signal, and consequently are not 
intended to imply (^timum sensor design criteria. The major results of this study can 
be summarized as follows. 

(•} A layered water balance model was developed which could 

be used to predict daily soil water contents for I'cm increments 
in the 0~30 cm depth soil zone within approximately 0.02 cmV 
cm'^. The model was successfully tested on four data sets. Additional 
testing appears warranted, however, primarily to determine whether 
the storage capacity values (estimated here to correspond to a soil 
tension of 0.09 bars) apply in conditions different from those under 
which the four data sets were collected. 
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(!i) Diurnal soil water conf’ent changes appear quite large near the 

3 3. 

soil surface (0.125 cm /cm in 0~0.5 cm layer if moisture 

3 3 

content is about 0,175 cm /cm ) but decrease rapidly with depth 

so that at the depth of 2"*3 cm, the diurnal range of moistures 

3 3 

becomes less than 0.03 cm /cm'^. It should be remembered that 
these values were based on measurements of only one soil in four 
different months (March, April, Septerdber, Octc^er). Although 
the exact values may differ from case to cose, it is apf orent that 
considerable near~surface diurnal soil moisture variation con 
occur. This should be taken into account when estimating soil 
moisture by real~time methods, particularly thc^e which respond 
primarily to surface and near’^urface soil moisture. 

(iii) In the cose of bare soii, large soil moisture changes were found to 
be confined largely to the top 30 cm of the soil, based on data 
collected in Eudora silt loam. This was primarily because the water 
lost by evaporation originated near the surface, and the precipitation 
amounts were on the average not sufficiently large(in spite of an 
abnormally wet season) to penetrate into the deeper layers. Thus it 
is concluded that under conditions of frequent precipitation, the 
error in computed vertical moisture distribution in a bare soil caused 
by assuming zero evaporation from depths greater than 30 cm would 
be sufficiently small or negligible under conditions similar to those 
of the experiment . 

(iv) If the effective precipitation during a certain period exceeds 
the storage capacity of the lower zone of the soil profile (depths 
greater than 30 cm in this study), it may be assumed that the lower 
zone contains only the amount equal to the storage capacity at the 
end of the period. The length of time needed for the excessive 
water to drain out will vary, however, depending on the hydraulic 
properties of the soil (and on ground water level if present)., 

(v) Moisture content near the surface, surface configuration, and 
surface roughness all affected surface and profile temperature of a 
bare soil in addition to the time of day. The variations were 
largest at the surface. Although all the temperature fluctuations 


154 


observed of rhe surface were propagated to greater depths, the 
differences encountered there were smaller than those at the 
surface; data generated by the soil temperature simulation model 
indicated that at depths greater than 30 cm, diurnal soil temperature 
variations can be expected to be around 5*^C or less . 

(vi) Microwave brightness temperature of a bare soil varied diurnal ly 
in response Ko changing moisture and temperature profiles. Large 
soil moisture changes, such as rapid evaporation loss following 
precipitation, can cause greater diurnal brighiness tempa'ature 
changes than relatively large soil temperoture fluctuations. 

(vii) A detailed analysis of the microwave brightness temperatures and 
associated moisture and temperature profiles revealed that even at 
the lowest frequency (lp4 GHz ), the brightness temperature was 
affected primarily by the air/soil discontinuity and much less so 

by the subsurface soil properties. That is, the vertical soil moisture 
and tenq>erature profile variations, although generated so os to 
represent relatively extreme conditions (for an agricultural area), 
were not large enough to substantially influence the brightness 
temperatures. For all frequency/moisture/temperature profile 
coR^inations investigated, the major dielectric discontinuity 
existed between air and the first soil layer which rendered the 
surface and near-surface moisture content to have a major effect 
on the brightness temperotures . Furthermore, a consideration of 
dielectric properties suggests that subsurface soil moisture gradients 
would hove to be in excess of those commonly found in a cultivated 
soil in order to produce reflections comparable to those at the soil 
surface. Similar results were arrived at when effective power 
reflection coefficients were computed for coherent radiation, such 
as transmitted by radar. 
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U Is importani' to remember fhaf i'he "surface layer dominaHon" 
was es^abl^shed using a model which neglected scattering. A cursory 
examination of the factws Involved suggested that the most effective 
scafterers In ibe soil would be dielectrically homogeneous soil clods 
separated by air spaces; the optimum clod size would depend on the 
wavelength. Since it Is not known to which extent such clods occur 
under field conditions, the conclusions derived on the basis of the 
brightness temperature model should not be considered final . 

(vili) By combining water balance, hydrological, and microwave 
remote sensing concepts, an algorithm was developed to 
monitor soil water storage In the profile as a function of time 
and location. Suitable primarily for an agricultural area, the 
algorithm could be used for the entire soil profile or for 
Individual zones within the profile. It would require the 
following Input: hydrological soil group, stcxoge capacity 
of the profile or of Individual zcmes, crop type, precipitation 
amounts and durations/neteoro logical variables needed for 
computing actual evaporation (evapotranspiration), and a rainfall 
pattern as determined by the remote sensors. 

(xi) Although active and passive microwave sensors have demonstrated 
their sensitivities to soil moisture variations for bare and vegetated 
fields, additional experimental work is needed before an optimum 
set os sensor{s) parameters con be defined such that the effects of 
soil surface roughness and vegetation cover are minimized. 
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APPENDIX A. 


SOIL TEMPER/^ JURE SIMULATION MODEL 


The purpose of the soil temperature model is to predict diurnal variations 
of surface as well as subsurface temperatures of a bare soil . As stated in Chapter 3, 
the prediction is accomplished by (!) computing the surface soil temperoture os a 
function of time, and (ii) using the surface temperature as a boundary condition to 
solve the one-dimensionol heat diffusion equation. These two steps are described 
in the following sections. 


A.l DETERMINATION OF SURFACE TEMPERATURE 


The magnitude of surface soil temperature at any time is governed by the 
energy exchange at and near the soil surface. It is therefore possible to estimate 
the surface temperature by determining magnitudes of individual terms in the energy 
balance equation which, neglecting horizontal radiation terms, can be written as 


where 


Therefore 


Rj^ = H + LE + G, 



= S - L 


a o 
= (1 - a) S. 


f 

/ 



S -L = H+LE + G. 
a o 


(33) 


(Individual symbols are defined at the end of this section). 
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The longwave radiaHon contribuHon fo net radiation consists of two 
parts, L and L. . If energy emitted by the surface L is considered a function of 
air temperature, then according to the Stefan-Boltzman law. 



The energy emitted by a cloudless atmosphere L. has been determined to follow the 
relationship (Sellers, 1965) 

L. = e (0.165 - 0.000769RH) . 

For a cloudy sky a factor (1 - k^n) approximates the reduction of outgoing radiation 
(Sellers, 1965), Furthermore, Sellers (1965) derived a correction factor to account 
for the difference between surface and air (screen height) temperatures in computing 
L^. Combining the above formulas, can be computed as follows: 

) 

I = 4eoT^ (T - T ) +e (1 - k,n) (0.165 - 0.000769RH). (34) 

o a ^ 5 a ^ 1 

Sensible heat H can be computed from the aerodynamic equation; 

pC^ k^u (T - T ) 


then 


H 


(cS) 


— TT~^ 


The basic equaHon for la^en^ heal- flux LE calculaMon is given by 


LE 


0.622 p L D (e^ - ) a 

p 


(36) 


As explained by Foster (1972), surface vapor pressure e^ can be eliminated from 
Equation (36) by employing a finite difference form of the Clausius"Clapeyron 
equation 


e " e 
s 


so 


0.622 Le (T - T ) 
sa ' s a 


(37) 


If is added to and subtracted from the left side of Equation (37) and the 
equation is written in terms of (e^ - e^)^ fhe result can be substituted in Equation 
(36) to yield 


LE = 


0.622 L D a , 0.622 L e (T - T ) 

I sa ' s a 

+ e 

R T ( R 


sa 


- e 


(3B) 
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Parameter a in Equation (38) accounts for the reduction of potential 
evaporation of the soil after moistuie content decreases below a certain level. 

It was assumed that actual evaporation AE equals potential evaporation PE if 
moisture content is above soil storage capacity, and that the ratio A^PE 
decreases when moisture content is below the storage capacity. Therefore 

a = 1 if mQ 2 § 0.95SCq2 
•"02 

= otherwise. (39) 

^^-02 

For a homogeneous semi“infinite soi5, the soil heat flux at the surface G(0#t) 
con be expressed in the form (Sellers, 1965') 


G(0,t)=AT 4 ^ 

o 

provided that the soil is heated in a periodic manner so that surface temperature 
follows Equation (41); 

T (t) = T + AT sin (lit, (41) 

*s s o 

The amplitude term AT can be removed from Equation (40) by usif^ 

Equation (41). Furthermore, Equation (40) can be mon'pulated to express G in 
terms of Cfg“T^). The resulting equation will be 


Before fhe individual ferins of fhe energy balance equafion are confined 
to yield a formula for computafion, Equafion (34), (35), (38), and (42) will be 
replaced by simpler expressions: 


'■o=S(^s-Ta»*'=2 
G = (T^ - T^) + C7. 

Then hhe surface temperarure can be derermined from Equarion (43): 


T 

s 


T 

a 


+ 


+ C3 + 


(43) 


where; 


= 4 €0- , 

e (1 - k°n ) ( 0.165 - 0.000769RH ), 

pC D 


C = • P 
^3 R T ' 


^ 4 = 


0.38688 L D a e 


a 


so 


^ 0.622LDa(e^^-e^) 

'-s TP 


C^= 0.7071 (XC( 0)^/^ ) , 


0 > 


T - T 

C7 = 0.7071 ( XC (1))^/^ " ^s ^ 


0 ) 


( Individual symbols used in ihis secHon indicate the Following variables; 



ratio of actual and potential evaporation, dimensionless; 

soil heat capacity, cal/cmV°Q 

specific heat of air, ca|/g/°C; 

transfer coefficient, ctn/min; 

soil surface water vapor pressure, mb; 

saturated water vapor pressure of the air, mb; 

actual air vopor pressure at the height z, mb; 

soil heat flux, ly/min; 

sensible heat flux, ly/mti., 

vonKarman constant, dimensionless; 

cloud type coefficient, dimensionless; 

latent heat of vaporization, cal/g; 

longwave radiation emitted by the soil, ly/min; 

latent heat flux, ly/min; 

counter longwove atmospheric radiation, ly/min; 
outgoing longwave radiation, ly/min; 

3 3 

average moisture content in the 0"2 cm layer, cm /cm ; 

cloud cover in tenths, dimensionless; 

atmcspheric pressure, mb; 

gas constant, R = 2870.4 mb cm^/g/®K; 

relative humidity, percent; 

net radiation, ly/ min; 

absorbed shortwave solar radiation, ly/min; 

3 3 

average storage capacity in the 0~2 cm layer, cm /cm ; 

incoming shortwave solar radiation, ly/min; 

air temperature, °K; 

surface soil temperature, ®K; 

wind speed, cm/min; 

height of meteorological measurements, cm; 

roughness length, cm; 

soil albedo, dimensionless; 

infrared soil emlssivity, dimensionless; 





AT^ = surface temperature amplitude, ®C; 

^ = apparent soil thermal conductivity, cal/cm/min/°C; 

P = air density, g/cm > 

O' = S tefan"Boltz man constant, a = 8.17 !y/min/®K ; 

w = angular frequency, radians/time increment.) 


A. 2 DETERMINATION OF SUBSURFACE TEMPERATURES 

Assuming horizontal uniformity and heat transfer by conduction only, the 
change in soil temperature is described by (Sellers, 1965) 


6T (zA) 
Si 





(44) 


If D'(z) is not constant with depth. Equation (44) cannot be solved 
analytically. A numerical solution is available, however. Equation (44) can 
be written in a finite difference form as follows (Richmeyer and Morton, 1967) 


Tf^ -t| _ D'(z)(Ti, -2T| +t|., ) 

At (Az)^ 


(45) 


where 

2 

D'(z)= soil thermal diffusivity, cm /min; 
i = depth increment sequential number; 

j = time increment sequential number; 

T = soil temperature, °C; 

At = time increment, min; 

Az = depth increment in the soil, cm. 

Equation (45) can be solved provided both the initial temperature os a function of 
depth and the temperotures at the upper (soil surfoce) and lower boundaries ore 
known os a function of time. If the lost Az increment is located sufficiently deep 
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I 


below fhe surface, fhe lower boundary condifion is constant. The upper boundary 
can be computed as stated in the previous section. The initio! temperature profile 
can be measured or orrived at by simulation; the latter procedure was described in 
section 3.2. 

A. 3 THE COMPUTER PROGRAM 

Figure 70 is a generalized flow chart of the computer program for soil 
temperature computations. The governing time parameters are PERL {= length 
of the period between measurements of meteorological parameters) and DT(=At). 

The number of intervals NHA for which sol! surface temperature T^ will be 
calculated is given by 

NHA = 1440/PERL. 

For the calculations described in this report, the values used were; PERL=10, 

DT=2.5, MHA=144. 

The soil variables needed are soil texture and moisture content, bulk density, 
field capacity, and initial temperature profile as a function of depth. The depth" 
dependent parameters are specified for NL layers; the layers were assumed to be 
2 cm thick in this study. 

The site is described by latitude, declination slope and aspect. Besides \ 

air temperature, relative humidity and wind speed at PERL increments, magnitudes > 

th . . 

of air pressure, surface roughness length, surface temperature for 0 time increment, ^ 

and of the average daily surface temperature are required. | 

The computation for a given time increment proceeds as follows: | 

— hwr angle; | 

— incoming solar radiation (see below); I 

— albedo (from data given by Idso et a| . (1975b)); | 

“ soil thermal conductivity and capacity (see below); | 

— surface temperature from Equation (43); | 

— profile temperatures in NTINC iterations using Equation (45);NTINC=PERl/DT; | 

—• evaporation from individual layers in the 0 to 30 cm zone, | 

based chi the method described in section 2. 2. 1.1 and latent I 

heat computation from Equation (38). i 
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The computoMon of surface temperature requires knowledge of soil thermal 
conductivity CONS, thermal capacity CAPS, and water content near the surface 
WCS (Equation 42, 43). These quantities were calculated as an overage of values 
at the surface (computed by assuming constant gradient for the 0 to 4 cm depth) 
and at the depth of 2 cm. 

Solar radiation incident on a generally sloping surface was computed as 
Q = QDR + QDF , 

where 


QDR = 


SMT (cos (ZEN) cos(SL) + sin(ZEN) sin(Sl) cos(AZS - ASP)) 

/ 

DIS^ ( 1 + (TC) (OM)) 


and 

QQp_ 0.38TC (cos (0»5Sk)) (SMT) 
DIS^ (1 + ( TC) (OM)) 


where 

SMT 

— 

solar meteorological constant, ly/min; 

DIS 


radius vector of the sum, dimensionless; 

TC 

= 

transparency coefficient for Kastrow’s equation 

OM 

— 

(Sivkov, 1971), dimensionless; 

actual optical mass of the atmosphere, dimensionless; 

ZEN 


zenith angle (0° = solar noon), degrees; 

SL 

= 

slope of the surface, degrees; 

AZS 


azimuth angle of the sun, degrees; 

ASP 

= 

aspect of the surface (0® = south-facing), degrees. 

Derivations of 

the formulas and values of the various coefficients used in the solar 


radiotion computations can be found in Sivkov (1971), Seilers (1965), List (1951), 
and Kondratyev (1969), 
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Thermal conducl-tvify was catculal'ed separately fc^ each layer using the 
procedure described by de Vries (1963; see also Wierenga et al ,, 1969). The 
basic formula is 


n 

Xrf i (^) 

X = 

where 

n = the number of different types of particles; 

X. = volumetric fraction of the i particles; 

^ = thermal conductivity of the particles. 

The subscript i = 0 refers to the continuous phase. Following de Vries (1963), the 

3 3 

continuous phase was water if moisture content was above 0.03 * IT cm/ cm (IT = 

1 fw sand, 2 for loam, 3 for clay) and air if moisture content of the layer was below 
this value. Besides air and water, two other kinds of particles were considered, 
namely quartz, and other mineral constituents; the reason for separating soil mineral 
fraction into two groups was the difference in X. between them. Values of A. for 
all components were taken from de Vries (1963) for the temperature of 20°C; 
temperature dependence of A. was not accounted for 'n the computations since it is 
generally small . The thermal conductivity of air-filled pores was equal to the sum 
of thermal conductivity of dry air, and of which accounts for the heat move- 
ment by vapor phase across the gas-filled pores (de Vries, 1963). The value of A^^ 
was cmstant, while \ was constant only above the field capacity; below field 
capacity, \ decreased linearly to A^ for a dry soil . The values of k. were 
computed as follows: 



k, = 4- 


(47) 


+ (X. - 


^ )g- 

O ' 

.th 


where g. is a factor that depends on the shape of the i particles, with g| + g 2 + 93 “ 
1 . For solid particles, the values used were ~ 92 ~ 83 “ 0*750 (Wierengd 

et ai., 1969). Values of g| and 22 for oir particles were assumed to decrease linearly 
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from 0.333 for wafer-saturofed soil fo some value of field capacify and fhen linearly 
to 0.013 for a dry soil; the value at field capacity depended on both porosity and 
field capacity in the given soil layer. For moisture contents below 0,03 * IT, the 
soil thermal conductivity was determined by interpolation between the \ value at 
0,03 * IT and X for a dry soil; in accordance with the data by Skaggs and Smith 
(1967), the X of a dry soil was first multiplied by a factor of 1,65, 

thermal capacity C was computed for each layer from the formula 


C = 0,2 Pg + m. 


(48) 


where 

pg = soil bulk density; 
m = moisture content. 


appendix b 


EXPERIMENTAL PROCEDURES 

Dafra from fhree different experiments were used in various ports of this 
study. Procedures employed to obtain these data are summarized in this Appendix. 


B.l EUDORA SOIL MOISTURE REGIME STUDY 


This experiment was carried out on an agricultural fielcT (latitude 38^57. 8*N, 
longitude 95°08.TW; NWNENE36-12S-20E) approximately 4 miles east from 
Lawrence, Kansas. The experiment was initiated following wheat harvest and plow- 
ing and progressed until tensiometric measurements had to be discontinued due to 
freezing temperatures. 

The experimental site was located on the floodplain of the Kansas River; 
consequently, the soils exhibit considerable spatial variability. Two soil mapping 
units were identified at the site, namely Eudora silt loam and Kimo silty clay loam 
(Dickey and Zimmerman, in press). The first soil occupied a somewhat elevated 
ridge, while the second one was located in a depression; no sharp boundary existed 
between them. Typical profiles of these soils have been described as follows (Dickey 
and Zimmerman, in press), 

B.T.I EudaroSllt Loam 

In a representative profile the surface layer is very dark grayish brown silt 
loam about 12 inches thick. The upper part of the substratum is dark grayish brown 
coarse silt loam about 36 inches thick. The lower part of the substratum is grayish 
brown coarse silt loam thinly stratified with a few sandy and clayey layers less than 
1 inch thick, 

’ The cooperation of Mr. A, Wichman is gratefully acknowledged. 
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Eudora soils are very friable and easily worked. They have high nafural 
ferhMify and available water capacity. Permeability is moderate. 

Nearly all the acreage of these soils is cultivated. Eudora soils are well 
suited to all crops commonly grown in the county. Also vegetable crops produce 
well on these soils. 

Representative profile of Eudora silt loam, 125 feet south and 50 feet west of 
the northeast corner of Sec. 15, T. 12 S ., R. 19 E ., in a cultivated field: 

Ap 0 to 7 inches, very dark grayish brown (lOYR 2 / 2 ) silt loam; 

moderate medium granular structure; very friable; mildly alkaline; 
gradual smooth boundary. 

A1 7 to 12 inches, very dark grayish brown (lOYR 3/2) silt loam; 

moderate medium granular structure; very friable; many worm 
casts; mildly alkaline; gradual smooth boundary. 

Cl 12 to 23 inches, dark grayish brown (lOYR 4/2); coarse silt loam; 
massive; very friable; many worm casts; mildly alkaline; gradual 
smooth boundary. 

C2 23 to 40 inches, dark grayish brown (lOYR 4/2); coarse silt loam; 

massive; very friable; few worm costs; at 33 to 34 inches a layer of 
silty clay loam; mildly alkaline; gradual smooth boundary. 

C3 40 to 43 inches, dark grayish brown (10Y8 4/2); silt loam; massive; 
very friable; mildly alkaline; clear smooth boundary. 

C4 48 to 72 inches, grayish brown (lOYR 5/2); coarse sil': loam; 

massive; very friable; sorrre thin, less than 1 inch thick of sandy 
and clayey layers; mildly oikaline; strong effervescence . 

The A horizon ranges from 10 to 20 inches in thickness. The A horizon 
may be dark grayish brown in areas where 6 inches or more of silty sediments were 
deposited by the flood in 1951. The A horizon can be coarse silt loam, loam or 
medium silt loam. The C horizon is usually a coarse silt loam with some layers more 
sandy or clayey textures. Below 36 inches the texture can be coarse silt loam, loam, 
very fine sandy loam, fine sandy loam or loamy very fine send that is usually 
calcareous. 


B.1.2 Kimo SUty Cioy Loom 

Represenl’Otive profile of Kimo stlfy clay loam, 250 foeh north and 60 feet 
west of southeast corner of sec. 7, T. 12 $., R, 20 E., in a cultivated field; 

Ap 0 to 6 inches, very dark gray (10 YR2|/1) heavy silty clay loam.! 
weak fine granular structure; firm; mildly alkaline; slight 
effervescence; gradual smooth boundary. 

A1 6 to 14 inches, black (10 YR 2/1) silty clay; weak fine and 

medium subangular blocky structure; very firm; mildly alkaline; 
gradual smooth boundary . 

A12 14 to 24 inches, very dark gray (10 YR 3/1) silty clay; few fine 

strong brown mottles in lower part; weak, fine, subangular blocky 
structure; very Hrm; mildly alkaline; clear smooth boundary. 

AC 24 to 28 inches, dory gray (10 Y? 4/1) light silty clay loam with 
some streaks and fingers of darker color; few, fine, strong“brown 
mottles; some weak fine granular structure and massive; friable; 
mildly alkaline; strong effervescence; gradual smooth boundary. 

lie 28 to 60 inches, grayish browm (10 YR 5/2) coarse silt loam; 

massive; very friable; miidly alkaline; strong effervescence. 

To monitor changes in soil moisture and meteorological parameters related 
to these changes, the following instrumentation was used: 

(i) Tensiometers were built using poly> h»yl chloride (PVC) pipe 
(Schedule 80, 1.27 cm I.D.), porous ceramic cups (manufactured 
by R&J Ceramics, 2552 Cottoge Way, Sacramento, California), 
acrylic tubing (1.59 cm O.D.), and translucent nylon tubing 
(0.107 cm I.D.). The porous cup and acrylic tubing were glued 
with contact cement to opposite ends of the PVC pipe cot to the 
required length. The nylon tubing was fit into a hole drilled 
through the PVC pipe wall and glued. During the meoairements, 
the upper end of the tensiometer was closed with a neoprene stopper. 
All tensiometers were tested for air leaks prior to installation. 

(ii) Rainfall gauge (diameter 10 cm). 

(iii) Hygrothermograph with a weekly recording cycle. 
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(iv) Anemomet'ers al- three heights (10 cm; 30 cm, 300 cm). 

(v) Pyranometer for solar radiation measurements; however, 
the recording equipment broke down permanently after two 
weeks of operation. 

The experiment was designed so as to encompass main sources of soil moisture 
variability, nomely depth, time, distance, and soil type. Depth was taken into 
account by installing the tensiometers 10 cm, 25 cm, 40 cm, 55 cm, 70 cm, 

85 cm, 100 cm, 115 cm, 130 cm, and 145 cm ("complete" set) or 115 cm, 130 cm, 
and 145 cm ("partial" set) below the soil surface. In either case, the nylon tubings 
from all tensiometem of one set v/ere fastened to a single meter stick and drew 
mercury from one container. One complete and one partial set were installed 2 m 
apart onihe Kimo silty clay loam (="Site 3"). On Eudora silt loam, two complete sets 
were installed 2 m apart ("Site 1”) and an additional partial set was placed approximately 
100 m away from the complete sets ("Site 2"). A dike (5 cm high above the soil 
surface) was installed around oil sets to prevent runoff. A neutron probe access tube 
was inserted at a 50 cm distance from each set; unfortunately, the neutron counts 
were found to give moisture contents Incompatible with both tensiometric and 
gravimetric measurements as well as the precipitation history, and therefore had to 
be discarded. 

All tensiometers were installed 30 and 31 July, 1973, Subsequently, mercury 
level reading and precipitation were taken at approximately 0730 every day (period 
S/2/73 to 9/6/73) or at intervals 1 to 3 days {9/6/73 to 11/3/73). 

Six (on Site 1) and three (on Site 3) undisturbed core samples 4 cm in diameter 
and 10 cm high were taken for moisture characteristics (i.e,, tension vs, moisture 
content relationship) determination. Due to unavailability of suction plotes,o 
Soil Moisture Equipment Corporation pressure membrane extractor was adjusted for 
measurements at low tensions by adding a water column monometer and an air pressure 
manometer (pressure range 0 to 60 p.s.1 .); water manometer was used fix* tensions up 
to 100 cm. Each set of somples was placed on a pressure plate (maximum pressure 3 
bars) into a larger basin; after slowly raising the water level within the basin so that 
the samples were almost submerged, the basin was covered to prevent evaporation 
from the samples, and the samples saturated for approximately 16 hours. Following 
this period, the pressure plate with the samples were transferred into a pressure 
cooker, closed tightly, and the pressure was set ot 0.01 atmosphere. Subsequently, 


air pressure was Increased to 0.02, 0.05, 0.1, 0,2, 0,5, 1 .0 and 2,0 atmospheres; 
when an equilibrium was established following every pressure increase, the outflow 
from the cooker was measured and recorded. At the end, samples were taken out 
weighed individually, dried ot 105^C, and weighed again. Moisture characteristics 
were then calculated from the finai moisture content, measured outflows, and from 
the pressures applied. Difficulties were experienced in processing samples from Site 
3, due mainly to their lower nurr^er and finer texture. Consequently, only moisture 
data from Site 1 are presented here. 

Tensiometric readings T (in cm) were converted Into values of suction S (in 
cm) using Equation 49*: 

S = 12.6 T" 13.6H - D + 14 , (49) 

where H = height of the mercury level in the container above soil surface in cm; 

D = depth of the porous cup below the surface in cm. 

Constant of 14 accounts for friction loss inside the nylon tubing. 

Using the moisture characteristic for a given depth, moisture contents 
corresponding to measured S values were determined for all measurements cm Site 1 
and 2. 


b.2 phoenix experiments 


Some data used in various parts of this study were generated during experiments 
carried out in 1970, 1971, and 1973 in Phoenix, Arizona, The experimental condi- 
tions we'e quite similar, however, and may therefore be summarized collectively. 

The following information was provided by Jackson (1973) and Idso et al , (1974). 

The enperiments were conried out on a 72 by 90 m field at the U, S. Water 
Conservation Laboratory in Phoenix, within which three lysimeters are installed. 

The soil, Adelanto loam (fihe-loamy, mixed, hyperthermic, Torrifluventic 
Haplustoll), is fairly uniform to about 100 cm and has been cultivated numerous times 


*Courtesy of Dr. W. Powers, Department of Agronomy, Kansas State University, 
Manhattan, Kansas. 


during the past years. At the start of each experiment, the lysimeters and surround” 
ing area were irrigated with approximately 10 cm of water. 

During the drying period, soil samples for gravimetric moisture content 

determination were taken at frequent intervals for depths of 0~0.5 cm, 0.5-1 cm, 

1-2 cm, 2-3 cm, 3-4 cm, 4-5 cm, 5-7 cm, and 7-9 cm; the sampling was not 

begun until some drying had occurred, thus allowing access to the field. The 

July, 1970, experiment was conducted between 10 July and 17 July, and sampling 

was done at 0.5 hour intervals. In the March, 1971, experiment {March 3 to 

April 8) measurements were taken from 5-18 March, 25 March, and 8 April, also 

at 0.5 hour intervals. In the 1973 experiments, samples were collected at 20 

minute intervals. Jackson et a| . (1973) estimated that the totol error in resulting 

3 3 

moisture content values was about + 0.001 cm /cm . Evaporation loss was measured 
by lysimeters at the same time os soil water content. 

Soil temperatures were measured by thermocouples at the soil surface and at 
several depths below the surface. For example, in the September, 1973 experiment, 
the interval between adjacent sensors was 0.5 cm (at depths less than 10 cm), 2 cm 
(10 to 35 cm), 15 cm (35 to 65 cm), and 35 cm (above 65 cm). 

In addition to soil temperature and moisture, several meteorologicol variobles 
were measured: incoming and reflected solar radiation, net radiation, air temperature, 
atmospheric vapor pressure, and wind speed . 


